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Abstract 


Experimental  investigations  and  theoretical  modeling  of  the  energy  deposition  and 
acceleration  processes  in  IkW-class  arcjet  thrusters  using  hydrazine  propellant  are 
presented.  The  report  highlights  detailed  research  in  arcjet  fundamentals  that  was 
reported  in  the  two  Ph.D.  theses  of  S.  A.  Bufton  and  T.  W.  Megli.h2 

We  report  for  the  first  time  the  evidence  of  thermal  non-equilibrium  (different 
electron  and  gas  temperatures).  Chemical  kinetics  for  hydrazine  dissociation  and 
ionization  play  important  roles  in  the  overall  energy  conversion.  Electrostatic  probe 
measurements  at  the  exit  plane  validate  the  model. 


^S.  A.  Bufton,  "Exit  Plane  Plasma  Measurements  of  a  Low-Power  Hydrazine  Arcjet," 

Ph.D.  Thesis,  Department  of  Mechanical  and  Industrial  Engineerins:,  Universitv  of 
Illinois,  1995.  6  6  > 

^T.  W.  Megli,  'A  Nonequilibrium  Plasmadynamics  Model  for  Nitrogen /Hydrogen 
Arcjets,"  Ph.D.  Thesis,  Department  of  Mechanical  and  Industrial  Engineering, 
University  of  Illinois,  1995. 


MEASURED  PLASMA  PROPERTIES  AT  THE 
EXIT  PLANE  OF  A 1  kW  ARCJET 


Abstract 

With  recent  advances  made  in  low  power  (1-2  kW)  arcjet  numerical  modeling, 
validation  of  numerical  models  through  experimentation  has  become  increasingfy 
important.  We  describe  multiple  electrostatic  probe  techniques  for  characterizing  flow 
conditions  at  the  exit  plane  of  a  1  kW-class  hydrazine  thruster.  An  improved  electrostatic 
(madruple  probe  tecimique,  which  is  also  modified  to  account  for  plasma  gradients  over 
the  finite  probe  volume  and  multiple  ion  species,  provides  measurements  of  ne,  Te  and 
ui/ Cm  at  the  thruster  exit.  Quadruple  probe  results  are  also  presented  for  measurements 
at  several  centerline  locations  immediately  downstream  of  the  thruster  exit  plane.  A 
spatially-resolved  time-of-flight  electrostatic  probe  technique  is  employed  for 
measurements  of  radial  profiles  of  the  plasma  axial  velocity  ui.  ITie  quadruple  probe 
results  (Te,  ne,  ui/ Cm)/  coupled  with  independent  measurements  of  plasma  velocity,  are 
iised  to  estimate  ion  temperature  Ti=Tg  and  the  extent  of  thermal  non-equilibrium  at  the 
thruster  exit.  Exit  plane  experimental  data  for  uj,  Te  and  Tg  are  compared  with 
cornputatiqnal  arcjet  model  predictions,  showing  agreement  for  ui  and  T^.  The  model 
underpredicts  Te  and  a  direct  comparison  of  ne  results  is  precluded  oy  large  axial 
gradients  at  the  thruster  exit. 
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Nomenclature 

Electrode  geometric  surface  area  [m^j 
Most  probable  thermal  speed  [m/s] 
Electron  charge  [coulombs] 

Ion  composition  parameter  (=nH+/ne) 
Probe  electrode  current  [A] 

Arqet  operating  current  [A} 

Ion  saturation  current  density  [A/m^] 
Boltzmann  constant  Q/K] 

Particle  mass  [kg] 

Propellant  flow  rate  [mg/s] 

Particle  density  [cm"3] 

Thruster  power  [kW] 

Designation  for  electrode  "k" 

Collision  cross  section  [m2] 

Electrode  radius  [mm] 

Electrode  clearance  [mm] 

Temperature  [K,  eV] 

Ion  flow  velocity  [m/s] 

Particle  thermal  speed  [m/s] 

Probe  electrode  voltage  [Volts] 

Electrode  relative  bias  voltage  [Volts] 
Floating  electrode  voltage  [Volts] 

Plasma  potential  [Volts] 

Arget  operating  voltage  [Volts] 

Axial  distance  from  thruster  exit  [mm] 
Wake  effect  collection  area  parameters 
Flow  divergence  angle  [Degrees] 
Dissociation  energy  [eV] 

Mathematical  Gamma  function 
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1.  Introduction  and  Background 

Although  research  in  the  area  of  electrothermal  thrusters  is  relatively  mature,  only 
recently  has  the  arcjet  become  applied  to  commercial  satellites.!  With  these  thrusters 
begmning  to  find  a  niche  in  industry,  accurate  numerical  models  become  increasingly 
importaint  as  design  tools  to  help  in  the  imderstanding  of  fundamental  arq'et  physics.  As 
such,  adequate  measurements  of  fundamental  plasma  properties  must  be  available  for 
the  purpose  of  validating  numerical  efforts.  It  is  the  purpose  of  this  research  to:  1) 
provide  fundamental  data  for  validation  of  H2/N2  computational  arcjet  models;  2)  add 
to  the  existing  base  of  fundamental  data  in  the  veiy  near-field  plume  region;  and  3) 
demonstrate  the  feasibility  of  electrostatic  probing  techniques  as  useful  diagnostics  in  the 
small  geometries  associated  with  the  1  kW  arcjet.  In  this  paper  we  present  electrostatic 
quadruple  probes  results  for  the  exit  plane  of  a  1  kW  hydrazine  arcjet. 
Where  applicable,  comparisons  are  made  with  experimental  data  and  computational 
model  predictions  performed  by  our  group. 

Electrostatic  probes  have  found  much  use  in  experimentally  determining  plasma 
parameters  in  space  and  laboratory  plasmas.  Several  authors^-^  have  employed  classical 
Langmuir^  single  probes  in  the  characterization  of  low  power  arcjet  plumes.  While  single 
probes  are  appealing  due  to  their  relatively  simple  construction,  interpreting  the  probe  V- 
^  difficult.  A  viable  alternative  is  the  triple  electrostatic  probe,^ 

which  allows  simultaneous  measurement  of  Te  and  ne  without  the  necessity  of  a  voltage 
sweep  and  a  corresponding  probe  V-I  characteristic.  The  utility  of  triple  probes  for 

has  been  demonstrated  in  MPD  thruster  plumes^'^^and  in  a  1 
kW  hydrogen  arcjet  plume.!3  Additionally,  a  general  review!^  of  electrostatic  probe 
techniques  and  their  uses  in  electric  propulsion  devices  has  recently  been  presented. 


Under  certain  conditions,  the  crossed-probe  technique!^"!^  can  be  used  to  measure 
(=plasma  velocity / ion  inost  probable  thermal  speed)  in  flowing  plasmas.  The 
quadruple  electrostatic  probe^^  combines  the  crossed^robe  and  triple  prooe  techniques 
for  ^multaneous  measurements  Te,  ne,  and  ui/ Cm.  Tne  utility  of  the  quadruple  probe 
has  been  demonstrated  in  an  MPD  thruster  and,  more  recently,  in  the  very  near-field 
plume  of  a  1  kW  hydrazine  arcjet.21 


11.  Quadruple  Probe  Technique 

Figures  la.  and  lb.  show  an  electrical  schematic  and  potential  plot  of  the  quadruple 
probe,  consisting  of  three  cylindrical  electrodes  (Pj,  P2,  P3  )  which  are  aligned  with  the 
^  j  flow  vector  and  one  electrode  (P4)  which  is  perpendicular  to  the  plasma  flow.  P3 
and  P4  are  biased  at  constant  voltages  Vd3  and  V(i4  relative  to  Pi,  and  P2  assumes  the 
floating  potential  Vf  of  the  local  plasma  (l2=0).  The  circuit  comprised  of  electrodes  1,  3 
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and  4  is  electrically  floating  such  that  ion  current  collected  at  P3  and  P4  is  balanced  by  the 
electron  curr^t  collected  at  Pi,  thereby  determining  the  potential  difference  Vd2 
l7Vd2(Te)J-  The  three  aligned  electrodes  1-3  are  similar  to  the  triple  probe  and  yield 
simultaneous  measurements  of  ng  and  Tg.  The  addition  of  the  peipendicular  electrode  P4 
allows  the  measurement  of  the  parameter  ui/cm,  determined  from  the  ratio  I4/I3  for 
crossed  electrostatic  electrodes  with  Vd3  =  Vd4-  For  the  quadruple  probe,  the  measured 
quantities  are  the  electrode  currents  I3  and  I4  and  the  potential  difference  between 
electrodes  1  and  2,  Vd2*  Currents  I3  and  I4  are  determined  by  measuring  voltage  drops 
across  resistors  (R3  =  R4  =  100  Q).  The  electrode  bias  voltages  Vd3  and  Vri4  are  provided 

by  batteries,  with  Vd3=Vd4=  12V. 


It  has  been  shown^  that  probe  contamination  can  have  a  marked  effect  on  Tg 
measurements  through  drastic  changes  in  the  measured  quantity  Vd2-  In  the  present 
study,  the  probe  was  cleaned  thoroughly  before  each  data  acquisition  session  using  ion 
bombardment  cleanmg.l3  xhis  was  accomplished  by  biasing  all  four  electrodes  24  V 
below  facility  ground  to  attract  ions  to  the  probe.  The  electrodes  were  then  placed  in  the 
arcjet  plasn^  stream  for  several  -Is  intervals,  so  that  the  electrodes  would  glow  and 
burn  off  me  contaminant  layer.  Quadruple  probe  sweeps  acquired  before  cleaning 

always  yielded  electron  temperatures  that  were  higher  than  the  post-cleaning  values  bv  a 
factor  of  ~  2.  r  o  j 

Probe  Construction 

A  schematic  of  the  quadruple  probe  used  in  this  study  is  shown  in  Fig.  2.  The  probe  is 
comprised  of  4  individual  tungsten  wire  electrodes  of  0.25  mm  diameter  with  an  exposed 
length  of  2.5  mm  for  electrodes  1-3  and  2.0  mm  for  electrode  4.  The  location  of  the 
perpendicular  electrode  has  been  modified  from  that  of  previous  probes, to  coincide 
with  the  axial  location  of  the  probe  center.  Each  wire  is  mounted  in  round  single-bore 
alumina  (AI2O3)  tubing,  which  is  in  turn  mounted  in,  and  supported  by,  tmn-wall 
stamless  steel  tubing.  This  assembly  is  mounted  in  a  larger  stainless  steel  tube  to 
tacihtate  connecting  the  probe  tip  to  the  probe  positioning  mechanism  discussed  below. 
Inside  this  larger  tubing  is  a  four-bore  alumina  tube  that  electrically  insulates  the  four 
probe  leads.  A  high  temperature  ceramic-based  adhesive22  is  used  to  bond  each  of  the 
probe  el^trodes  and  components  in  place.  A  small  quantity  of  adhesive  is  placed  at  the 
junction  between  P4  and  its  single  bore  alumina  tube  to  electrically  insulate  part  of  the 
electrode,  thereby  reducing  the  error  associated  with  determining  its  geometric  surface 
area  (Fig.  2.).  °  ® 

The  c^ter-to-center  separation  between  electrodes  1  and  2  is  1.4  mm.  This  spacing, 
along  with  electrode  lengths  of  2.5  mm,  allows  an  unaligned  quadruple  probe  to  handle 
flow  divergence  angles  6  ~  30  without  "shadowing"  of  one  electrode  by  the  other.  In  an 
effort  to  minimize  the  effect  of  the  probe  geometry  on  the  plasma  flow,  a  gradual 
trarisition  between  the  electrodes  and  the  probe  support  has  been  implemented.  The 
radius  of  the  single-bore  alumina  electrode  supports  is  ~  0.4  mm,  which  minimizes 

?erturbations  of  the  plasma  flow  at  the  junction  between  the  electrodes  and  their  support. 

he  exj^sed  length  of  the  single-bore  alumina  is  ~  5  mm.  Ceramic  adhesive  is  usecfto  fill 
the  void,  and  ease  the  transition,  between  the  electrodes  and  the  stainless  steel  probe 
support  tubing.  ^ 

Probe  Length  Scales 

For  the  conditions  expected  at  the  exit  plane  probe  location  (Te~0.6  eV,  ne  ~  4  x  lO^^ 

''  Debye  length  is  2.9  x  10"^  cm  and  the  ratio  of  electrode  radius  to  Debye  length 
rp/A,D=45,  so  that  the  ion  sheath  surrotmding  the  electrodes  can  be  classified  as  thin  with 
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probe  measurements  5  mm  downstream,  this  ratio  increases  to  rn/^D 
-  55.  The  clearance  (s)  between  electrodes  1  and  3  (Fig.  2)  is  s~l  mm  so  that  s/Xu  >  250 
and  the  sheaths  on  adjacent  electrodes  do  not  interact. 

To  avoid  strong  sensitivity  of  the  ion  current  to  small  misalignments  between  the  flow 
^  axis,  the  so-called  end  effect  parameter23 

l“^*"/ui^D)(kTe/mi)-i/^)]  should  be  greater  than  50.  For  our  plasma  conditions,  this 
parameter  is  conservatively  estimated  at  Tl>250.  The  electrocles  should  also  be  long 
enough  that  the  Bohm  sheath  ion  current  collected  along  the  cylinder  length  (Ib)  is  much 
larger  than  the  current  collected  at  the  probe  end  due  to  the  convection  of  charged 
particles  into  the  probe  tip  (It).  A  first  orcier  analysis  of  the  ratio  Ib/It  for  the  conditions 
at  the  exit  plane  showed  Ib/It  >  30,  so  that  the  error  introduced  to  the  collected  ion 
current  by  neglecting  the  cylinder  tip  is  less  than  4%. 

Ihie  radius  rp  of  the  tungsten  electrodes  (Fig.  2)  was  chosen  as  a  compromise  between 
the  thm  sheath  assumption  (rp/^D  »1)  and  the  assumption  of  free  molecular  flow  over 
the  probe  electrodes  (Amfp/rp  »1).  For  the  seven  species  plasma  expected  at  the  exit 
f  11  '  relevant  collision  mean  free  paths  were  calculated, 

he  followmg  assumptions  are  made  in  calculating  mean  free  paths:  Te=7000  K,  Tg=2500 
X  10  cm"  ,  and  nneutral  =  3  x  10^^  cm^  with  10%  dissociation  at  the  thruster 
following  (worst  case)  mean  free  paths  were  calculated:  Xn-n/rp  >20,  Aj-n/rp 
>5,  Ax-i/rp  >30,  Xi-e/rp  =20,  Xe-n/ rp  >20,  and  Xe-e/ rp  =30.  Clearly,  the  probe  electrodes 
are  operating  in  the  collisionless  regime.  Additionally,  Xmfp/XD  >200  for  all  collisions, 
with  an  overwhelming  majority  of  the  collisions  having  Xnxfp/Xo  >1000,  indicating  that 
the  collisioniess  sheath  assumption  is  appropriate. 

Quadruple  Probe  Theory 

Quadruple  probe  analysis  assumes:  1)  The  ion  sheath  surrounding  the  tungsten 
wires  IS  thin  compared  to  the  wire  radius  (rp  »  Xd);  2)  Each  probe  wire  and  sheatii  are 
collisioniess  (Xixxfp  »  rp  »  Xd);  3)  The  electron  energy  distribution  is  Maxwellian;  4) 
Velocity  slip  between  plasma  species  is  negligible;  and  5)  A  wake  region  of  low  ion  and 
electron  density  exists  downstream  of  the  perpendicular  electrode  P4. 

Derivation  of  quadruple  probe  response  is  similar  to  that  of  the  triple  probe^  and  has 
been  presented  elsewhere.  In  our  previous  work21  we  presented  some  of  the 
revisions  necessary  to  allow  the  application  of  triple  and  quaaruple  probe  theory  to 
(N-H)  plasmas.  In  the  present  study,  the  prone  is  applied  to  a  simulated 
hydrazine  plume,  with  multiple  species  of  ions  present,  whereas  most  of  the  previous 
applications  of  triple  and  quadruple  probes  have  been  for  plasmas  with  one  crominant 
species  of  positive  ion.  Single  Langmuir  probes  have  been  used  in  nitrogen  and 
arcjet  plumesA"^  but  their  use  is  not  contingent  on  knowledge  of  the 

In  the  present  study,  we  present  a  significant  modification  of  the  quadruple  probe 
meory.  Rather  than  utilizing  an  analytic  expression  for  the  current  collected  at  electrode 
T  T  current  collected  at  P3  through  the  ratio  of  the  measured  quantities 

I3  and  I4.  The  significance  of  this  modification  is  tnat  it  completely  removes  errors  in  Te 
and  rie  measurements  that  were  due  to  assumptions  regarding  effective  current  collection 
arM  of  P4.  For  a  quadruple  probe  with  electrode  geometric  surface  areas  A3  (=Ai=A2) 
and  A4  and  electrode  biasing  such  that  Vd3=Vd4,  the  following  expres-sion  can  be 
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(1) 


1  _  1  +  exp((|)Vd3)  -  2exp[(|)(Vd3-Vd2)  ] 

{I4/I3 )  {exp[(t>(Vd3-Vd2)  ]  -  1 } 

Since  Vds  is  prescribed  and  Vd2  is  experimentally  determined,  Eq.  (1)  determines  Tp 
=e/k(j)  through  iteration. 

The  ion  saturation  current  density  ji  is  related  to  measured  quantities  by 

(I3/A3)  (I  +  I4/I3) 
exp((j)Vd2)-l 

^  (2)  is  used  to  determine  ne  from  I3, 14  and  Vd2  by  relating  ji  to  the  electron  density 
through  the  Bohm  sheath  analysis.  For  a  plasma  with  more  than  one  species  of  collected 
ion,  the  Bohm  analysis  must  be  modified  to  account  for  the  collection  of  each  ion  species 
Bohm  sheath.  For  the  region  of  the  hydrazine  plume  of  interest,  H+  and  N+  are 
the  dominant  ions  species,  and  the  Bohm  expression  becomes:21 


ji  =  e(kTe)l/2exp(-i)(_JlH^ 


nN+  \ 

(mN+)^^^  I 


(3) 


For  ne=nH‘''  +  nN+,  Eqs.(2)  and  (3)  yield 

,  (4) 

e(kTe)^/2  [  exp(<t)Vd2)  -  1  ] 

Eq^  (4),  K  represents  the  effect  of  multiple  ion  species  on  the  probe  response  to  electron 


where  |i  =  (mH+/ mN+)l/2.  For  a  plume  contain-ing  only  H+,  nN+  =  0.0  and  K  reduces  to 
umty,  such  that  Eq.  (4)  reduces  to  single  component  (H+)  form. 

Experimental  observations24/25  show  that  low  nitrogen  ion  densities  are  found  near 
the  exit  plane  in  aminonia  and  hydrazine  arcjets.  It  is  suggested  that  this  is  a  result  of 
charge  exchange  collisions  in  the  expanding  nitrogen-hydrogen  propellant.  It  can  be 
shown  (Fig.  3)  that  equilibrium  number  density  ratios  PN/nH  and  nN+Zne"*"  are  equal 
<  7500  K  at  p  =  1  atm.  For  this  temperature  range,  nisj/nH<  0.4,  due  to  the 
pr^erential  dissociation  of  H2  (£d  =  4.48  eV)  over  N2  (ed  =  9.61  eV).  TTie  arcjet  exit  plane 
N-H  plasma  is  highly  dissociated  and  partially  ionized,  with  Te  >  Tg,  and  the  estimated 
r^ge  of  gas  temperature  is  Tg  <  6000  K.  Ionization  by  electron  collision  is  negligible,  and 
nitro^n  rapidly  undergoes  charge-exchange  collisions  with  hydrogen,  the  cross-section 
which  IS  on  the  same  order  of  magnitude  as  N-H  momentum  transfer  collisions, 
inus:  H  +  N"*"  < — >  H+  +  N,  and  in  the  absence  of  other  ionization  processes,  the  two 
charge-exchange  reactions  balance  in  steady  state.  Writing  these  volumetric  collision 
rates  in  the  form  nQv,  with  relative  thermal  speeds  vH  N  ~  vH,  and  cross  sections  QH- 
N+  =  QN-H+  gives: 


»N+-nN 
nH+  nH 


(6) 
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which  directly  implies  a  low  nitrogen  ion  density  in  the  exit  plane  (Fig.  3).  Thus,  the 
nitrogen  ion  density  is  effectively  determined  by  the  gas  temperature  Te  and  not  the 
electron  temperature.  For  Tg  <  6000  K,  nN+/nH+  <  0.1  and  K  <  1.07,  indicating  a  ~7% 
mcrease  of  ne  due  to  the  presence  of  nitrogen  ion  probe  current. 

Crossed  Electrostatic  Probes 

The  use  of  crossed  electrostatic  probes  for  measuring  plasma  velocityf^'f 7  ^nd  the 
implementation  of  the  crossed  probe  technique  into  the  quadruple  probe  have  been 
present^  elsewhere,  por  a  cylindrical  electrode  oriented  at  an  angfe  0  relative  to  the 
/  Yf /f®  collected  ion  current  is  a  function  of  the  quantity  (ui/cin)sin0,  where  Cm 
-  (2kTi/mi)-i/^.  For  the  case  of  two  equally  biased  probes,  with  one  probe  aligned  with 
the  flow  yector  and  the  other  one  normal  to  it,  the  single  species  collected  ion  current 
ratio  for  the  thin  sheaths  case  is:  17 


I4  _  2  A4 
I3  V5fA3 


]  2 

n=0 


L  n!  J 


r(n+2.) 

2 


(7) 


^  order  to  apply  the  crossed  probe  technique  to  the  multicomponent  (H2-N2)  plasmas  of 
interest  here,  significant  revisions  must  be  made  to  Eq.  (7).  In  the  original  derivation  of 
the  crossed  probe  technique  a  Maxwellian  velocity  distribution  is  assumed  for  the  single 
component  ion  collection.  15  For  the  present  study,  this  expression  has  been  rederived 
using  a  two-component  form  of  the  Maxwellian  velocity  distribution.  The  derivation  is 
similar  to  the  single  component  case,  and  results  in  a  contribution  to  the  ratio  I4/I3  from 
each  of  the  collected  ion  species.  For  the  two  species  of  ions  assumed  present,  the  ion 
composition  is  defined  by  f  =  nH+/ne  and  (1-f)  =  nisj+/ne  (ne=ni-i+  +  The  resulting 
expression  for  the  measured  current  ratio  is 


k 

13 


n=0 


+M.(xN*)(l-f)exp 


where  the  constant  C  is  given  by 


(8) 


2A4 

V5f  A3[f +  (l-f  Jp] 


(9) 


^1^*4-  ^  '  the  first  term  in  the  brackets  { }  is  the  H+  contribution  and  the  second  term  is 
me  INI  contribution.  The  most  probable  thermal  speed  for  N"*"  (cni,N+)  h^s  been  written 
in  t^ms  of  C]:n,H+  (since  Cin,N+=  P-  Cm,H+)-  Note  that  each  of  the  terms  is  weighted  by  the 
relative  density  of  its  corresponding  ion  [f,  (1-f)] . 

Because  a  wake  is  formed  behind  P4,  a  portion  of  the  electrode  area  A4  does  not 
collect  10ns.  Tne  parameters  xh+  and  xn+  are  the  fractions  of  the  electrode  area  A4 
available  mr  H  and  N+  collection,  respectively.  Because  the  extent  of  the  wake  effect  is 
determmed  by  the  relative  magnitudes  of  the  directed  ion  velocity  uj  and  thermal  speed 
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of  the  ions,  xh+  and  xn+  are  not  necessarily  equal.  For  Tg=Ti=2500  K  and  a  directed 
velocity  ui  ~  6500  m/ s,  uj/ Cin,H+  ~  1  and  ui/ Cm,N+  ~  4.  Several  references  discuss  the 
wake  effect  for  high  speed  flows.l^'i^  It  has  been  postulated  and  experimentally  verified 
for  an  argon  plasma  at  uf/ Cm  ~  1,  that  the  wake  effect  causes  the  effective  collection  area 
to  be  half  of  the  geometric  electrode  area.l^  Since  ui/Cin,H+  is  expected  to  be  ~  1,  xh+  is 
assumed  to  be  1/2  so  that  only  the  front  half  of  P4  collects  ions.  For  the  nitrogen  ions,  the 
wake  effect  reduction  in  collection  area  is  expected  to  be  more  pronounced.  Since  the 
directed  velocity  ui  is  several  times  Cin,N+/  N+  is  assumed  to  be  collected  by  the  projected 
area-*^"  of  the  perpendicular  electrode  only,  such  that  XN+=l/7r. 

Given  f,  Eq.  (8)  becomes  an  expression  for  ui/cin,H+  versus  the  measured  current 
ratio  I4/I3.  Equation  (8)  is  plotted  in  Fig.  4  for  0  <  f  <  1,'xh+=1/2  and  xN+=l/7r.  Pure  H+ 
collection  is  represented  by  f=1.0,  2H+  +  N+  collection  by  f=0.67,  and  pure  N+  collection 
by  f=0.0. 

If  an  independent  measure  of  the  directed  plasma  velocity  Ui  is  known,  F^.  4  can  be 
used  to  determine  Cin,H+  [=(2kTi/mH+)l/2]  and,  hence,  the  gas  temperature  T^=Ti.  We 
describe  below  uj  measurements  which  lead  to  an  estimate  or  Ti. 

Effect  of  Radial  Gradients  on  Probe  Response 

Quadruple  probe  measurements  at  the  thruster  exit  plane  show  smooth  symmetric 
profiles  for  ne  and  the  electrode  currents  I3  and  I4,  as  snown  below.  However,  Vci2  is 
known  to  have  a  very  as^^metric  profile,^!  which  is  interpreted  as  a  result  of  the  steep 
radial  gradient  in  ne  off-axis,  and  to  a  lesser  extent,  the  radial  gradient  in  Te-  These 
gradients  are  perpendicular  to  the  quadruple  probe  electrodes  and  the  arcjet  thrust  axis. 
In  the  presence  of  plasma  radial  gradients,  the  individual  electrode  current  equations^O 
are  rewritten  in  terms  of  the  plasma  parameters  at  each  electrode,  and  are  rearranged  to 
solve  explicitly  for  Vd2  and  I3  in  terms  of  radial  profiles  of  Te  and  ne: 


V2=V2  (Te2,1C) 

(10) 

Vi=Vi  (Tel,  Te3,  riel,  ne3, 14/I3,  K) 

(11) 

Vd2=V2-Vi 

(12) 

13=13  (Te3/  ne3/  K) 

(13) 

The  subscripts  on  the  parameters  Te  and  rie  in  Eqs.  (10)-(13)  denote  the  local  value  of 
those  parameters  at  each  of  the  electrodes  (1-3)  of  the  quadruple  probe.  The  ratio  I4/I3  is 
a  mea^red  quantity  and  K  is  a  function  of  the  plasma  composition  discussed  above. 
Smce  the  centerline  experimental  data  are  at  zero  gradient  conditions,  these  data  can  be 
used  with  Eqs.  (1)  and  (4)  to  determine  the  centerline  values  of  Te  and  ne- 

Using  ^s.  (10)-(13),  Fig.  5  displays  Vd2  vs.  radial  position  for  the  quadruple  probe  of 
yig-  2,  for  fiat  and  Gaussian  ne  and  Te  profiles.  For  Gaussian  ne  and  constant  Te,  Vd2 
mcreases  monotonically  across  the  thruster.  For  a  Gaussian  Te  distribution  broader  than 
ne,  the  Vd2  profile  shows  more  asymmetry.  The  constant  Te  case  can  be  solved 
analytically  for  Vd2,  giving: 


In  this  equation  the  effect  of  the  density  gradient  can  be  clearly  seen.  As  the  swept  probe 
starts  outside  the  plume  and  enters  it,  ne3/ riei  <  1.  At  the  axis  nei  =  ne3,  and  as  the  probe 
leaves  the  plume,  ne3/nei  >  1.  The  result  [Eq.  (14)]  is  a  monotonically  increasing  Vd2 
across  the  thruster  face. 
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Given  ^itrary  profiles  of  Te(r)  and  ne(r),  we  see  that  Vd2(r)  and  l3(r)  can  be 
predicted.  The  inverse  problem,  of  extracting  ng  and  Tg  profiles  from  measured  Vd2  and 
13,  can  also  be  solved.  Since  yd2  [Eq.  (12)]  is  a  strong  function  of  Tg,  and  I3  [Eq.  (13)]  is  a 
strong  function  of  ng,  there  is  sufficient  uncoupling  that  an  iterative  routine  is  used  to 
£-1  ^  ^e(f)  and  ng(r)  profiles  whicn  reproduce  the  measured  Vd2  and  I3 

profiles.  Starting  with  centerline  values  and  assumed  Gaussian  profiles,  the  routine 
marches  the  probe  position  across  the  thruster  face,  while  generating  a  new  Tg  profile. 
This  updated  profile  is  used  with  Eq.  (13)  to  update  ng(r),  and  the  process  is  repeated 
until  both  profiles  converge.  Typically  five  numerical  iterations  are  required  to  reach 
convergence. 


III.  Time-of-Flight  Velocimetry 

Since  in  addition  to  the  quantities  ne  and  Te  a  quadruple  probe  measurement  yields 
the  quantity  ui/ (2kTi/mH+)^/,  knowledge  of  the  heavy  particle  velocity  ui  can  be  used 
to  determine  the  hea\w  particle  temperature  Tg  =  ij.  W^e  have  designed  and 
implemented  a  time-qf-flight  velocirnetry  technique  for  measuring  spatially  resolved 
axiai  velocities  in  arcjet  plumes.  This  technique  is  derived  from  a  previous  Current 
Velocimetry  (CMV)  method. ^3,26  in  our  method,  a  short  duration  current 
deficit  is  superimposed  on  the  arcjet  current,  the  effect  of  which  is  monitored  by  an 
electrostatic  time-of-flight  (T-O-F)  probe  as  it  convects  axially  in  the  thruster  plume. 

Time-of-Flight  Probe 

T^e  CMV  technique  can  be  made  spatially  resolved  with  electrostatic  probes,  thereby 
eliminating  the  need  to  Abel-invert  the  data  as  required  by  line-of-sight  optical  methods. 
Figure  6  shows  the  electrostatic  time-of-flight  sensmg  probe  used  in  tnis  study.  It  consists 
of  two  0.75  mm  diameter  tungsten  wires  which  are  electrically  insulated  along  most  of 
their  length  by  two  1.6  mm  o.d.  alumina  tubes.  The  exposed  length  of  each  tungsten 
electrode  is  0.25  mm  and  the  axial  separation  between  the  two  is  5  mm  (±  2%).  The 
elevation  of  the  probe  is  such  that  the  arcjet  thrust  axis  is  located  between  the  sensing 
electrodes.  The  electrodes  are  biased  24  V  below  facility  ground  (arqet  anode)  to  monitor 
the  local  electron  saturation  current  density  rie,  each  measurement  location. 

Variations  in  probe  electrode  saturation  currents  are  monitored  as  they  convect  over  the 
two  T-O-F  electrodes,  with  the  temporal  separation  of  these  signals,  combined  with  the 
probe  geometry,  indicating  the  flow  velocity. 

Although  this  probe  is  physically  intrusive,  the  probe  electrodes  (0,75  mm  dia.)  are  in 
me  near-free  molecular  flow  regime  and  are  thus  expected  to  have  little  effect  on  the 
flow.  Additionally,  any  plasma  perturbations  due  to  the  electrodes  (shocks,  collisional 
effects,  etc.)  are  expected  to  have  a  similar  result  at  each  electrode,  thereby  having 
minimal  effect  on  the  convective  time  delay. 


Arcjet  Current  Deficit  Circuitry 

power  processing  unit  (PPU)  produces  current  with  a  high  frequency  ripple 
^  15-20  kHz,  ±  8%  about  the  mean)  which  causes  'natural"  fluctuations  in  ne  anci  Te- 
These  inherent  fluctuations  are  difficult  to  utilize  with  the  T-O-F  probe  due  to  their  long 
period  (--  55  |is),  compared  with  the  expected  convective  time  delay  of  the  T-O-F  probe 
1  \is).  Instead,  fluctuations  are  artificially  introduced  to  the  arcjet  plasma  flow  by 

superimposing  a  short  duration  (--  a  few  |is)  current  deficit  pulse  on  the  arqet  operating 
current.  ;  r  o 


8 


T^e  circuit  used  to  modify  the  arc  current  is  similar  to  that  used  in  previous  CMV 
studies  and  is  shown  in  Fig.  7.  The  current  pulse  is  initiated  by  closing  a  switch  to  a 
fast  R-C  circuit  in  parallel  with  the  arcjet,  causing  a  sudden  decrease  in  arc  current  while 
the  capacitor  charges.  This  short  duration  (~  4  )xs)  current  deficit  results  in  a  "tagging"  of 
the  plasma  in  the  arc-heating  region.  Switch  closure  is  accomplished  with  an  IRF-SSO 

u  transistor  (FET),  which  is  closed  by  the  gate  pulse  shown  in  Fig.  7.  Originally 

the  R-C  circuit  element  was  replaced  by  a  single  resistor  (20-30  Q),  but  ringing  on  the 
edges  of  the  resulting  square  current  pulse  obscured  the  desired  T-O-F  probe  signals. 

Spatially  resolved  axial  velocities  are  measured  by  triggering  the  arcjet  current  deficit 
jmlse  at  several  different  T-O-F  probe  locations.  Since  the  probe  is  swept  radially 
through  the  plume  during  the  experiment,  it  is  necessary  to  trigger  the  current  deficit 
pulse  repeatably  at  each  desired  radial  location.  Probe  location  is  measured  by 
monitoring  the  voltage  of  a  precision  potentiometer  (Fig.  7)  which  has  been  calibrated 
against  the  arcjet  anode  for  probe  position.  The  position  (voltage)  of  the  desired 
measurement  location  is  set  at  the  voltage  comparator  circuit,  so  that  a  trigger  signal  is 
issued  to  the  HP  214A  pulse  generator  when  the  probe  reaches  the  desired  radial 
location.  In  turn,  the  pulse  generator  supplies  a  square  pulse  to  "close"  the  FET  and 
initiate  the  arcjet  current  deficit.  Because  1)  the  gate  pulse  must  be  positive  (relative  to 
the  FET  voltage),  and  2)  the  arcjet  cathode  (and  FET)  operates  at  — 100  V  (anode 
gnunded),  the  FET  Driver  Circuit  (Fig.  7)  is  a  necessary  buffer  between  the  FET  and  the 
HP  fiulse  generator.  Since  the  probe  is  swept  relatively  slowly  (~  20  cm/s),  the  probe 
location  is  effectively  "frozen"  during  the  ~  3-4  |lis  elapsed  time  between  the  initial  trigger 
and  the  monitoring  of  the  tagged  plasma  at  the  T-O-F  electrodes. 

IV.  Experimental  Apparatus 

Electrostatic  probe  experiments  were  carried  out  at  the  exit  plane  of  a  1  kW  NASA 
arcjet  thruster  operating  at  a  flow  rate  of  50  mg/s  of  2H2  +  N2  to  simulate  fully 
decomposed  hydrazine.  Each  component  of  the  propellant  was  individually  metered  by 
Unit  rnass  flow  controllers  and  mixed  in  the  propellant  line.  During  steady  state 
operation,  the  thruster  operated  at  10.0  A  and  112  V  (p/rh  =  22.4  MJ/kg).  Arc  current  was 
measured  with  a  Hall  effect  current  transducer.  The  thruster  exhausted  into  aim 
diameter  x  1.5  m  long  vacuum  tank^l  at  ~200  mTorr  during  arcjet  operation.  Anode 
temperatures  were  monitored  with  an  optical  pyrometer,  reaching  ~  1200  K  during 
steady  state  operation.  ^  ry  o  & 

The  quadmple  probe  and  T-O-F  probe  are  accurately  positioned  and  repeatably  swept 
throum  the  mruster  plume  with  the  probe  moimt  system  shown  in  Fig.  8.  Tlie  complete 
a^embly  is  fastened  to  a  step  motor  (Fig.  8)  mounted  on  a  linear  translation  carriage.^1 
The  motor  and  probe  assembly  carriage  is  capable  of  ±  5  cm  linear  translation 
perpendicular  to  the  ardet  axis,  allowing  the  probes  to  be  manually  swept  through  the 
plurne  at  ~  20  cm/s.  Carriage  and  probe  position  are  determined  to  0.25  mm  by  a 
precisiori  10  kD  linear  potentiometer  coupled  to  the  carriage  manual  drive  shaft,  and 
calibrated  against  the  outer  diameter  of  the  anode.  The  probe  support  arm  is  coupled 
direcuy  to  the  motor  shaft  with  the  probe  tip  on  the  shaft  axis,  allowing  probe  angle  to  be 
varied  (m  0.9  increments)  without  changing  the  probe  tip  location. 

Probe  elev^ion  is  repeatably  located  to  within  0.5  mm  by  aligning  the  probe  center 
with  a  line  etched  on  the  arqet  anode.  The  separation  between  the  probe  tips  and  the 
thruster  exit  plane  is  reliably  set  to  within  0.1  mm  with  a  spark  gap  gauge.  The  ardet 
mount  IS  known  to  deflect  toward  the  probe  ~0.7  mm  while  under  vacuum.  Any  error 
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associated  with  the  arcjet  exit  plane  position  because  of  this  is  systematic  and  does  not 
affect  the  relative  spacing  of  subsequent  axial  measurement  positions. 

V.  Experimental  Results  and  Analysis 

Experimental  quadruple  and  T-O-F  probe  results  are  presented  for  the  very  near-field 
plume  of  the  arcjet.  To  quantify  the  extent  of  axial  gradients  in  Tg  and  ng,  quadruple 
probe  measurements  are  made  at  several  probe  tip  locations  from  1-5  mm  from  the  exit 
plane.  Radial  profiles  with  probe  tips  located  1  mm  from  the  exit  are  presented  for  both 
the  quadruple  and  T-O-F  probes. 

Quadruple  Probe  Results 

For  the  quadruple  probe,  Tg  and  rig  are  determined  by  measurements  of  Yd2, 13  and 
14.  Since  these  parameters  are  expected  to  be  influenced  by  the  inherent  current  ripple 
associated  with  the  arcjet  PPU,  a  discussion  of  the  effect  of  the  ripple  on  the  probe  is 
warranted.  Figure  9  shows  the  centerline  ion  current  density  js  (=  13 /A3)  measured  at 
e^ctrode  3  ana  the  arc  current  larc  versus  time.  The  current  measured  at  electrode  3  is  in 
phase  with  the  arc  current,  except  for  a  delay  of  ~  3-4  ps  associated  with  the  time  required 
for  the  plasma  to  flow  from  the  arc-heating  region  near  the  cathode  to  the  quadruple 
probe.  Accounting  for  differences  in  average  velocity  and  measurement  location,  tnis 
delay  is  consistent  with  previous  triple  probe  results  for  the  hydrogen  arcjet.l3 

relationship  between  centerline  measurements  of  Vd2  and  j3  is  shown  in  Fig.  10. 
While  j3  closely  follows  the  phase  of  the  PPU  current  ripple,  Vci2  is  180°  out  of  phase  with 
P-  Smce  j3  ~  ng(Tg)l/2  and  Vd2  increases  monotonically  with  Tg  [Eq.  (1)],  Figs.  9  and  10 
indicate  that,  except  for  the  convective  time  delay,  ng  increases  and  Tg  decreases  with 
increases  in  larc  during  a  typical  PPU  ripple  cycle.  Specifically,  Eqs.  (1)  and  (4)  indicate 
that  the  ±  10%  current  ripple  causes  an  associated  ±  18%  ripple  in  ng  (in  phase)  and  a  ± 

5  /o  ripple  in  Tg  (180  out  of  phase)  about  their  means.  For  the  quadruple  probe  centered 
from  the  exit  plane,  the  mean  centerline  values  are  Tg  -0.6  eV  and  ng  -3.6  x 
10  cm"  .  The  Tg  results,  although  out  of  phase  with  larc/  are  in  phase  with  the  arc 
voltage  Varc  (which  has  a  ±  3%  ripple)  due  to  the  negative  impedance  characteristic  of  the 
arqet.  Thus  the  Tg  behavior  is  consistent  with  the  variation  of  the  E-field,  which  heats  the 
electrons.  These  results  are  different  from  previous  hydrogen  arcjet  data,f3  which 
indicate  that  ng  decreases  and  Tg  increases  with  increasing  Igro  It  is  unclear  at  this  time  if 
this  IS  an  indication  of  a  fundamental  difference  in  the  operating  physics  of  the  hydrazine 
and  hydrogen  arcjets. 

Because  of  the  large  variation  in  the  measured  quantities  Vd2/ 13  and  I4  with  the  PPU 
^rrent,  the  mean  values  of  these  parameters  are  used  to  determine  the  mean  values  of 
Tg,  ng  and  ui/  Cni,H+  through  Eqs.  (1),  (4)  and  (8).  For  radial  profile  data,  a  moving 
average  smoothing  routine  is  applied  to  the  raw  data  so  that  the  effects  of  the  PPU  ripple 
are  eliminated. 

Tg  and  ng  Axial  Profiles 

Previous  no  data  indicate  that  axial  gradients  over  the  length  of  the  quadruple  probe 
may  be  large.4.28,29  jo  quantify  the  effects  of  gradients  in  Tg  and  ng  over  the  2.5  mm 
quadruple  probe  length,  centerline  data  were  acquired  for  several  probe  axial  locations. 
Five  probe  tip  locations  from  1  to  5  mm  from  the  thruster  exit,  in  1  mm  increments,  were 
utilized.  Results  for  the  measured  current  density  j3  for  five  centerline  axial  locations  are 
shown  in  Fig.  1 1 .  The  geometry  of  the  three  aligned  electrodes  is  also  shown,  with  its 
placement  corresponding  to  the  probe  location  for  the  first  axial  location.  A  50  % 
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decrease  in  13  is  indicated  over  the  5  mm  region  investigated,  which  represents  the 
combmed  effect  of  changes  in  Tg  and  ng,  since  j3  ~  ng(Tg)l  /2 . 

Figure  11,  along  with  the  correspondii^  Vd2  data  and  Eqs.  (1)  and  (4),  yields  the  axial 
centerline  profiles  of  Te  and  n^  shown  in  fig.  12.  As  a  check  of  the  quadruple  probe  and 
to  insure  that  the  presence  of  the  perpendicular  electrode  does  not  influence  the  Tg  and 
ng  data,  independent  measurements  were  made  with  a  triple  probe  at  the  first  axial 
location  (x=  2.25  mm).  The  quadruple  and  triple  probe  Tg  and  rig  data  were  consistent 
within  5  %,  well  within  the  experimental  error  associated  with  the  experiment.  This  is 
not  a  surprising  result,  since  in  the  limit  I4/I3  ->  0,  Eqs.  (1)  and  (4)  reduce  to  the  triple 
probe  equations  for  Tg  and  ng. 

The  electron  temperature  (Fig.  12)  varies  from  6700  to  3600  K  (-0.6  -  0.3  eV)  between 
2.25<  X  <6.25  mm  from  the  thruster  exit.  These  results  are  roughly  consistent  with 
pevious  far-field  plume  results  of  Tg  ~0.2-0.3  eV  for  30<x<150  mm  and  Tg  -  0.1-0.2  eV 
for  x=320  mm.'^'^  Further,  the  present  data  seem  to  corroborate  previous  estimates^S  of 
Tg  >  0.5  eV  in  the  arcjet  nozzle  based  on  arguments  that  Tg  is  on  the  order  of  the  atomic 
excitation  temperature. 

Results  for  the  electron  density  axial  profile  are  also  shown  in  Fig.  12.  Over  the  5  mm 
range  of  probe  locations  ng  droj^  by  -  30  %,  with  the  magnitude  of  dug/dx  increasing 
slightly  near  the  thruster  exit.  The  ng  results  of  the  present  study,  along  with  those  of 
other  researchers,  are  summarized  in  Fig.  13.  While  the  ng  gradient  based  on  the  present 
data  is  large  near  the  exit  (dng/ dx  —  0.5  x  10^2  cm"^/mm  at  x=  2-3  mm),  it  is  small 
compared  with  the  axial  gradient  indicated  in  the  nozzle  interior  and  at  the  exit  plane.  It 
thus  appears  likely  that  large  axial  variations  in  the  ng  gradient  (large  d2ne/dx2)  are 
aresent  m  the  thruster  very  near-field  plume.  The  present  ng  data  appear  consistent  with 
previous  results  with  the  exception  of  previous  data  for  x=0  and  x=10  mm  (Fig.  13^29 
Towever,  these  data  represent  upper  limits  placed  on  rig  based  on  Stark  broadening  of 
me-of-sight  (non- Abel  inverted)  emission  spectroscopy  signals,  and  thus  are  consistent 
with  our  probe  results. 

Based  on  the  axial  variations  of  Tg  and  rig  reported  above,  it  is  apparent  that  the 
^adruple  probe  results  represent  plasma  conditions  averaged  over  the  axial  dimension 
(2.5  nim)  of  the  probe.  Additionally,  the  very  large  ng  gradients  presumed  to  exist 
iminediately  downstream  of  the  thruster  exit  require  extra  care  in  locating  and  reporting 
probe  positions  for  both  optical  and  intrusive  measurement  techniques. 

Exit  Plane  Radial  Profiles 

Traces  for  the  measured  quadnmle  probe  voltage  Vci2  and  current  I3  are  shown  for  a 
typical  exit  plane  radial  sweep  in  Fig.  14.  These  data  have  been  smoothed  so  that  the 
traces  shown  represent  the  mean  signals  (with  the  PPU  ripple  effects  removed).  As 
discussed  above,  the  I3  trace  is  symmetric  about  the  thruster  centerline,  while  the  Vh? 
trace  e^ibits  the  asymmetry  attributed  to  radial  gradients  in  ng  and,  to  a  lesser  extent,  Tg 
over  the  probe  face.2l  The  modified  quadruple  probe  theory  [Eqs.  (10)-(13)]  is  used  to 
solve  for  the  radial  profiles  Tg(r)  and  ng(r)  resulting  from  the  radial  quadruple  probe  data 
Vd2  and  I3  (Fig.  14).  This  analysis  indicates  that  the  Tg  profile  is  very  flat  near  the 
^  nim)  and  nearly  inform  at  Tg  -  0.6  eV.  This  uniform  Tg,  along  with  Eq. 
(13),  IS  used  to  determiiie  the  radial  ng  profile  at  x=2.25  mm  from  the  thruster  exit,  shown 
in  Fig.  15.  Uncertainty  in  rig  values  far  from  the  centerline  is  larger  due  to  misalignment 
between  electrode  3  and  the  flow  divergence  angle  8. 

T-O-F  Velocity  Probe  Results 
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The  effect  of  the  time-of-flight  current  pulse  on  thruster  operation  is  demonstrated  in 
Fig.  16,  which  shows  the  arqet  current  larc  and  voltage  Varc  versus  time  during  a  t5mical 
pulse.  The  current  drops  suddenly  as  the  FET  switch  closes  (Fig.  7)  and  the  0.1  |iF 
capacitor  begins  to  charge.  As  the  capacitor  nears  the  end  of  the  charge  cycle  (~  a  few  1 
P-s  R-C  time  constants),  the  current  through  the  FET  leg  of  the  circuit  decreases  so  that  the 
arc  current  resumes  the  normal  PPU  ripple  after  a  small  positive  ring.  The  width  of  the 
arqet  current  deficit  pulse  is  ~  4  |xs.  The  nature  of  the  behavior  of  Varc  is  not  immediately 
clear  from  Fig.  16.  However,  the  derivative  of  the  arc  current  (dlarc/ dt)  has  the  same 
qualitative  behavior  as  Varo  implying  a  possible  inductive  effect  resulting  from  the  arqet 
geometry  and/ or  from  the  current  pulse  circuitry. 

The  T-O-F  probe  response  is  shown  in  Fig.  17  for  a  typical  centerline  axial  velocity 
measurement.  The  signal  measured  at  the  downstream  electrode  has  been  increased  by  a 
factor  of  6  to  account  for  the  lower  electron  current  (lower  Te  and  rie)  5  mm  downstream 
of  the  upstream  electrode.  Both  the  upstream  and  downstream  electrodes  show  noise 
beginning  at  t  ~  2  |xs  that  corresponds  with  the  action  of  closing  the  FET  switch. 
Approximately  3  |is  later,  the  effect  of  the  arcjet  current  pulse  is  monitored  as  a  decrease 
in  the  electron  current  measured  by  the  upstream  electrode.  After  a  time  delay  for  the 
tagged  flow  to  traverse  the  5  mm  electrode  separation,  a  similar  decrease  is  noted  at  the 
downstream  electrode.  For  both  of  the  signals,  the  width  of  the  decrease  in  the  measured 
electron  current  is  ~  4  fxs,  approximately  the  same  as  that  of  the  current  deficit  pulse. 

Axial  velocities  are  inferred  from  the  convective  time  delay  demonstrated  in  Fig.  17 
and  the  known  separation  (5.0  mm)  between  the  two  T-O-F  electrodes.  The  time  delays 
are  derived  from  the  probe  response  using  the  temporal  separation  between  the  minima 
(due  to  the  decrease  in  electron  current)  measured  at  each  electrode.  This  is 
accomplished  by  taking  the  derivative  of  each  (smoothed)  signal  as  shown  in  Fig.  18.  The 
difference  in  the  locations  of  the  zero  crossings  is  the  time  required  for  the  tagged  flow  to 
travel  the  5.0  mm  electrode  separation. 

Results  of  T-O-F  probe  centerline  velocity  measurements  are  shown  in  Fig.  19  for  data 
acquire  during  three  individual  periods  after  thruster  ignition  (cold  start).  ITie  data 
u  centerline  axial  velocity  initially  averages  ~  6.0  km/s  during  the  period 

that  the  thruster  is  warming  up.  After  approximately  5-7  minutes  of  run  time,  the  axial 
averages  ~  6.5  km/s.  These  data  are  consistent  with  our  previous  results^l 
which  indicate  that  exit  plane  ne  measurements  reach  steady  state  after  ~  4-8  minutes. 
Durmg  the  warm-up  period,  the  measured  velocities  fluctuate  between  5.8  and  6.3  km/s. 
Smiilariy,  the  steady-state  centerline  axial  velocities  range  from  6.3  to  6.7  km/s.  TTiese 
variations,  although  approximately  within  the  experimental  error,  may  be  physical. 
Previous  work  utilizing  CMV  for  a  1  kW  hydrogen  arcjet^^  has  shown  similar  velocity 
fluctuations.  In  that  work,  it  was  determined  that  the  phase  of  the  PPU  ripple  was  not 
r^ponsible.  Similarly,  the  arcjet  current  deficit  pulse  was  always  triggered  at  the  same 
phase  in  the  PPU  currerit  ripple  (Fig.  16)  for  the  data  presented  in  this  study.  It  is  imclear 
what  may  be  causing  this  phenomenon. 

The  T-O-F  probe  was  also  used  for  off-centerline  measurements  to  generate  an  axial 
velocity  profile  at  the  thruster  exit.  For  these  data,  an  attempt  was  made  to  align  the 
probe  with  the  flow  divergence  angle  6.  It  was  found  that  using  a  fixed  angle  (0=0°) 
probe  produced  a  flat  or  even  inverted  profile  with  a  minima  on  the  centerline.  An  off 
probe  that  is  grossly  misaligned  with  the  flow  divergence  angle  (0«5)  will  see 
^  faster  plasma  at  the  downstream  electrode  than  at  the  upstream  electrode. 
Additionally,  the  electrode  separation  along  the  streamline  is  less  than  the  geometric 
separation  (5.0  mm)  by  the  factor  cos(5).  Both  of  these  factors  artificially  increase  the 
measured  axial  velocities,  with  the  former  phenomenon  having  a  greater  effect  than  the 
latter.  Although  the  plasma  streamlines  are  not  known  a  priori,  the  error  associated  with 
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T-O-F  probe  misalignment  for  off  axis  measurements  can  be  minimized  by  rotating  the 
probe  to  align  it  with  the  expected  flow  divergence  at  a  given  location.  The  flow 
divergence  has  been  shown"*  to  be  quite  large  ~  10  mm  from  the  thruster  exit,  and  is 
expected  to  be  large  at  the  exit  as  well.  For  T-O-F  measurements  at  r  =  0, 1,  2,  3  and  4 
mm,  the  probe  angles  used  were  0  ,  10.8  ,  18.0°,  25.2°  and  45.(f,  respectively.  It  was 
assumed  that  velocities  measured  in  this  manner  were  along  the  probe  angle  and  were 
converted  to  axial  components. 

Figure  20  shows  the  results  of  axial  velocity  measurements  made  at  the  thruster  exit 
plane,  with  the  upstream  electrode  located  1.0  mm  from  the  exit.  The  profile  is  fairly 
symmetric  and  is  compared  with  a  computational  model  as  discussed  below. 

Analysis  of  Crossed  Probe  Data 

Since  a  measurement  of  the  centerline  value  of  the  plasma  velocity  uj  has  been  made, 
quadruple  probe  data  can  be  used  to  estimate  the  heavy  particle  kinetic  temperature 
Tg=Ti.  For  the  Quadruple  probe  located  at  x=2.25  mm,  the  measured  ratio  of  the  ion 
currents  collected  by  electrodes  3  and  4  is  I4/I3  =  0.55  to  0.60.  This  current  ratio,  along 
with  Fig.  4  for  f=0.9  (K=1.07),  yields  ui/cm,H+  =  0.95  to  1.05,  which  indicates  that  the 
most  probable  thermal  speed  for  H+  (cm,H+)  is  approximately  equal  to  the  directed  ion 
velocity  ui.  For  the  measured  centerline  velocity  ui=6.5  km/s,  the  heavy  particle 
temperature  [Ti=(cmJH+FmH+/2]  is  Tg=Ti  ~  2500  K  (0.2  eV).  Therefore,  with  the 
previous  results  for  Te  ~  6700  K,  the  extent  of  thermal  non-equilibrium  on  axis  at  the 
thruster  exit  is  Tg/Tg  ~  0.4. 

Since  the  perpendicular  electrode  P4  in  reality  measures  I4  based  on  a  distribution  of 
Ui  and  Tf  along  its  length,  a  discussion  of  the  errors  introduced  by  the  probe  size  is  in 
order.  Because  the  quadruple  probe  theory  for  determining  Tg  and  rig  uses  the  measured 
value  of  J4/I3,  no  error  due  to  the  presence  of  the  perpendicular  electrode  is  involved  in 
determining  Tg  and  rig.  This  has  been  verified  (^ove)  with  independent  triple  probe 
measurements.  In  interpreting  the  crossed  probe  data,  however,  the  geometry  of  the 
probe  becomes  a  factor.  The  current  I4  is  likely  slightly  underrepresented  due  to  changes 
i^/cr®\i79  electrode  length.  However,  because  ui/cm,H+  is  proportional  to 

s  Tj  profile  peaked  at  the  centerline  has  the  opposite  effect  on  14.  Assuming  a 

4/  l4/l3=  1.1(0.57)~0.63.  The  ratio  ui/cm,H+  becomes  ~  1.15,  and  Ti  becomes 
2UUU  K  (0.17  eV). 

Comparison  of  Computational  Model  and  Data 

A  computational  arcjet  modeling  effort^O-31  been  undertaken  by  our  group 
concurrent  with  the  exit  plane  probe  studies  discussed  above.  TTie  model  is  an 
axisynimetric,  seven  species  hydrogen/nitrogen  plasma  code  that  utilizes  a  PISO 
algorithm  to  solve  the  computational  domain  up  to  the  exit  plane.  Separate  energy 
equations  are  formulated  and  solved  for  the  electrons  and  the  heavy  species.  The  anode 
tenmerature  distribution  is  included,  and  the  plasma  electrical  conductivity  is  coupled 
with  the  plasma  properties,  which  allows  a  self-consistent  solution  for  the  current 
distribution.  The  model  is  capable  of  both  chemical  equilibrium  and  non-equilibrium 
el^^  h  ^^^30  3^  complete  discussion  of  the  model  and  assumptions  has  been  presented 

Comparisons  are  made  between  experimental  results  and  model  predictions  for  a  1 
kW-class  arcjet  operating  on  50  mg/s  of  2H2+N2  to  simulate  fully-decomposed 
nydrazine  propellant.  The  arc  current  is  10.0  A  for  both  model  and  experiment.  The 
specific  powers  for  the  laboratory  thruster  (112  V,  10  A)  and  numerical  model  (94  V,  lOA) 
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are  22.4  and  18.2  MJ/kg,  respectively.  The  model  specific  power  is  contingent  on  the 
empirical  value  of  the  anode  fall  potential  used.  For  the  model  results  discussed  below,  a 
chemical  non-equilibrium  simulation  is  employed. 

pie  centerline  exit  plane  rie  predicted  by  the  model  is  ne=3  x  10l4  cm'^,  which,  given 
the  largely  varying  axial  ne  gradients  discussed  above,  seems  consistent  with  previous 
experimeiUal  data  (Fig.  13).  Because  of  these  large  gradients,  a  direct  comparison 
between  the  quadruple  probe  exit  plane  data  (probe  tip  =  1  mm  from  exit)  and  the 
computational  model  is  difficult. 

Centerline  exit  plane  heavy  particle  temperature  predictions  agree  favorablv  with 
experiment  (Tg^model  =  2800  K,  Tg^exp.  ~  2000-2500  K),  but  Te  results  differ  by  a  factor  of 
~  2  (Te,model  -  3000  K,  Te,exp.  ~  o6(J0  K).  Based  on  experimental  results  of  the  present 
and  previous  studies,  it  thus  appears  that  the  numerical  model  may  be  underpredicting 
electron  temperatures  at  the  thruster  exit  plane.  Note  that  model  predictions  or  thrust,  ui, 
Isp  and  Tg  are  unaffected,  smce  the  relative  electron  concentration  is  low  (ng/nn  <  0.001). 
Internal  nozzle  diagnostics^^  performed  by  our  group  are  expected  to  help  resolve  the 
discrepancy  between  experiment  and  numerical  model  predictions. 

^  comparison  between  the  predicted  axial  velocity  profile  and  the 
results  of  the  T-O-I^robe  surveys.  Agreement  is  excellent  and  within  experimental  error 
in  most  locations.  The  numerically  predicted  profile  is  slightly  more  narrow  and  peaked 
than  the  measured  profile.  This  is  an  expected  result  because  the  model  flow  is 
constrained  in  the  riozzle  (no-slip  at  wall),  which  differs  from  the  free  expansion  that 
occurs  m  the  plume  in  the  vicinity  of  the  T-O-F  probe. 

VI.  Summary  and  Conclusions 

]^perimental  results  for  two  electrostatic  probe  diagnostic  techniques  are  presented 
for  the  very  near-field  plume  of  a  1  kW  hydrazine  arcjet.  An  improved  quadruple  probe 
theory,  showing  excellent  agreement  with  independent  triple  probe  surveys,  was  used  to 
axial  profiles  of  Te  and  rie  for  x=2.25-6.2B  mm  from  the  exit  plane. 
While  the  gradient  in  ng  is  slight  oyer  the  region  investigated,  indications  are  that  the 
c^nsity  varies  by  1-2  orders  of  magnitude  within  the  first  few  millimeters  downstream  of 
the  exit  plane.  For  the  axial  position  nearest  the  thruster  exit,  Te=6600  K,  ne=3.6  x  10^2 
cm  and  Tg=Ti  ~  2000-2500  K.  With  the  heavy  particle  temperature  measurement,  the 
extent  of  thennal  non-equilibrium  has  been  quantified  at  the  thruster  exit.  A  spatially 
resolved  techni^e  for  measuring  plasma  axial  velocities  has  also  been  developed  and 
implemented,  '^e  radial  velocity  profile  was  measured  at  the  exit  plane,  with  a  peak 
centerlme  velocity  of  6500  m/ s.  The  above  results  were  used  to  evaluate  the  predictions 
of  a  numerical  H2/N2  arciet  model.  Agreement  was  very  favorable  for  both  the  exit 
j^ane  centerline  Tg  and  radial  uj  profile.  Based  on  previous  and  current  data,  it  appears 
®  may  underpredict  exit  plane  Te.  Internal  nozzle  diagnostics  are  expected 

to  illuminate  the  causes  of  this  discrepancy.  Comparison  of  ne  predictions  and 
experimental  data  could  not  be  attempted  due  to  the  large  gradients  in  rie  present  at  the 
exit  plane.  Clearly,  there  are  significant  advances  to  be  made  in  understanding 
fimdamental  arcjet  physics,  through  both  experimentation  and  modeling  efforts. 
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Fig.  la.  Quadruple  probe  electrical  schematic.  Measured  quantities  are  Vd2/ 13  and  I4. 


P2  P3  Pi  P4 

Fig.  lb.  Quadruple  probe  electrode  potential  plot  for  Vc}3=Vd4. 


Side  View  End  View 


Fig.  2.  Schematic  of  the  quadruple  probe  tip  swept  through  the  arqet  plume. 
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Fig.  3.  Ionic  and  atomic  number  density  ratios  for  simulated  N2H4  at  p=l  atm. 


f  (=  nH+/ne)  ~  0.9, 


Fig.  5.  Calculated  quadruple  probe  response  to  axisymmetric  Gaussian  profiles  of  Te  and 
ne  with  standard  deviation  o  (mm). 
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Fig.  6.  Side  view  schematic  of  the  time-of-flight  electrostatic  probe. 
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Fig.  8.  Schematic  of  the  probe  mount.  Probe  axes  can  be  rotated  in  0.9°  increments. 
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Time,  fis 

Fig.  9.  Arcjet  current  and  ion  saturation  current  density  showing  the  effects  of  the  PPU 
ripple.  The  delay  between  Igrc  and  js  is  ~  3-4  |is. 


Time, 

Fig.  10.  Vd2  and  j3  vs.  time.  Figs.  9  and  10  show  that  ne  increases  and  Tp  decreases  with 
mcreases  in  larc- 


Distance  from  Thruster  Exit,  mm 

Fig.  11.  Axial  variation  of  the  ion  saturation  current  density.  Horizontal  error  bars  are  2.5 
mm  long,  indicating  the  length  of  electrode  3. 
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Distance  from  Thruster  Exit,  mm 


12.  Axial  variation  of  Te  and  n^.  The  n^  gradient  is  -0.5  x  10l2  cnr^/mm  nearest 
the  thruster.  Experimental  error  in  Tg  is  ±  15%. 


Distance  From  Nozzle  Exit,  mm 

Fig.  13.  Comparison  of  centerline  ng  results  from  the  present  work  and  previous  studies. 


Radius,  mm 


Fig.  14.  Smoothed  profiles  of  l3(r)  and  VH9(r)  measured  by  the  quadruple  probe. 
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Radius,  mm 

Fig.  15.  Measured  rie  radial  profile  for  x=2.25  mm  from  the  thruster  exit  plane. 


Time, 

Fig.  16.  Variation  of  larc  and  Varc  during  the  current  deficit  pulse.  The  pulse  width  is 
much  smaller  that  the  period  of  the  PPU  ripple. 
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Time, 


^^^i^tion  of  the  T-O-F  probe  signals  during  the  current  deficit  pulse.  Probe 
electrode  signal  delay  is  ~  0.77  ps.  o  &  t 
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Time,  jis 

Fig.  18.  Derivatives  of  the  smoothed  T-O-F  probe  signals  depicting  zero-crossings  at 
signal  minima  locations.  i:  o  o 


Arcjet  Run  Time,  minutes 

ip*  Measured  centerline  axial  velocity  vs.  time  for  arcjet  warm-up.  The  numerical 
model  prediction  is  also  shown. 


6  4  2  0  2  4  6 


Radius,  mm 

Fig.  20.  Experimental  measurements  and  model  predictions  of  the  axial  velocity  profile  at 
the  thruster  exit. 
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A  PLASMADYNAMICS  MODEL  FOR  NONEQUILIBRIUM 
PROCESSES  IN  N2/H2  ARCJETS 


Abstract 


A  general,  two-dimensional,  steady,  non-LTE  (local  thermodynamic  equilibrium) 
moder IS  presented  for  nitrogen /hydrogen  arcjet  thrusters.  Both  thermal  and  chemical 
nonequilibrium  processes  are  incorporated  to  describe  a  multi-temperature,  multi- 
species  N2/H2  plasma.  The  current  and  anode  temperature  distributions  are  fully 
coupled  to  the  plasma  flowfield  with  realistic  boundary  conditions.  Thermal 
nonequilibrium,  chemical  equilibrium  results  are  presented  for  simulated  hydrazine. 
Itesults  for  both  thermal  and  chemical  nonequilibrium  hydrazine  and  hydrogen  arciets 
show  significant  thermal  nonequilibrium  near  the  anode,  and  high  frozen  flnw  losses. 
It  is  concluded  that  nonequilibrium  chemistry  models  are  necessary  to  capture 
accurately  the  interior  distributions  of  the  plasma  flowfield. 

Nomenclature 


9 

St 

go 

I 

Isp 


L 


kc 

Le 

Lh 

mi 

mij 

m 

M 

ni 


P, 


jj 

9ii 

Je 

Jg 


u 

V 

V 

V 
w 

x,r,0 

xi 

yi 


mean  thermal  speed  (m/s) 
coulomb  charge 

effective  diffusion  coefficientim^/s) 
elastic  energy  transfer  (W/m^) 

Mavitational  acceleration  at  sea  level  (m/s^) 
anck's  constant  (J-s) 
enthalpy  of  species  i(J/kg) 
heavy  species  total  particle  energy  (J) 
integrated  total  current  (A) 
specific  impulse  (s) 
current  density  (A/m^) 

Boltzmann's  constant  (J/K) 

Coulomb  constant  (Nm^/C^) 

Forward  rate  coefficient  (m^/s  or  m^/s) 
Reverse  rate  coefficient  (my/s  or  m^/s) 
Electron  inelastic  loss  rate  (W/m^) 

Heavy  species  inelastic  loss  rate  (W/m^) 
mass  of  species  i  (kg) 
reduced  mass  (kg) 
mass  flow  rate  (kg/s) 

Mach  number 


number  density  of  ith  species  (m"3) 
pressure  (N/ni2) 

mtegrated  ohmic  power  deposition  (W) 
collision  cross  section  (m2) 
anode  temperature  (K) 
electron  temperature  (K) 
h^v^  species  temperature  (K) 

axial  velocity  (m/s) 
radial  velocity  (m/s) 
velocity  vector  (m/s) 
voltage  (volts) 

azimuthal  swirl  velocity  (m/s) 
cylindrical  coordinates 
mole  fraction  of  species  i 
mass  fraction  of  species  i 
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^con  constrictor  length  (mm) 

Xmol  mixture  parameter  in  XinolN2  +  H2 
Zi  thermqd5mamic  partition  functions 
cci  catalytic  efficiency  of  species  i 
Ed  dissociation  energy  (J) 

£i  ionization  energy  (J) 

e  anode  emissivity 

anode  thermal  conductivity  (W/m-K) 

Ae  electron  thermal  conductivity  (W/m-K) 

Ah  heavy  species  thermal  conductivity  (W  /m-K) 
Ae  Spitzer  logarithm  term 

|i  gas  viscosity  (kg/m-s) 

V  collision  frequency  (s"l) 

6vH2  H2  vibrational  temperature  (K) 

6vN2  N2  vibrational  temperature  (K) 
p  mass  density  (kg/my) 

o  electrical  conductivity  (mhos  /  m) 

T.  stress  tensor  (N / m^l 

Hij  collision  inte^al  (m^) 


I.  Introduction 

Elechothermal  arcjets  offer  a  significant  advantage  in  specific  impulse  and  cost  over 
cony enhonal  satellite  propulsion  systems.  Low  power  1-2  kW  systems  have  been  flight 

J  •  used  for  north-south  station  keeping,  while  higher  power  20-30 

^  advanced  development  for  orbit  transfer  applications.^ 

The  arcjet  is  shown  schematically  in  Fig.  1 .  The  fluid  dynamics  and  energy  transfer  within 
the  arcjet  thruster  are  nonlinear  and  strongly  coupled.  The  propellant  is  injected  upstream  of  the 
constrictor,  typically  with  an  azimuthal  (or  'swirl')  velocity  component.  The  energy  conversion 
mechanism  for  an  arcjet  is  an  arc  discharge,  produced  by  a  voltage  difference  between  the 
mruster  nozzle,  which  is  the  anode,  and  a  conical  cathode  on  the  upstream  side  of  the  constrictor, 
typical  operating  voltages  for  1  kW  class  arcjets  are  100  V  DC  at  a  current  of  I  =  10  A.  The 
^c  current  distribution  depends  on  several  factors  in  addition  to  the  geometry.  The  distribution 
IS  coupled  to  the  propellant  mass  flow  rate,  composition,  thermal  properties,  electrical 
conductivity,  and  gas-dynamic  properties. 

Both  chemical  and  thermal  processes  in  the  arcjet  are  described  as  nonequilibrium.  Arc 
current  is  converted  to  electron  thermal  energy  through  ohmic  dissipation.  The  electrons 
transier  thermal  energy  to  heavy  species  in  the  arc  plasma  through  collisions.  This  energy  is  then 
converted  to  thrust  as  the  fluid  accelerates  tlyough  the  nozzle.  In  regions  of  low  pressure  and/or 
^  u  heating,  collisional  coupling  between  electrons  and  heavy  species  may  not 

est^lish  equal  gas  and  electron  temperatures.  Additionally,  flow  velocities  are  large  so  that 
iluid  residence  times  in  the  nozzle  are  of  the  order  of  a  few  ps;  much  of  the  energy  invested  in 
aissociation  and  ionization  is  frozen.  Pressures  drop  from  roughly  1  atm  at  the  constrictor  to 
miliiton:  conditions  at  the  exit;  radial  diffusion  of  electrons  from  the  arc  core  is  important  in 
determining  the  arc  structure. 

The  thrust  and  specific  impulse  produced  by  an  arcjet  are  determined  by  several  factors 
including  the  power  transferred  to  the  propellant,  the  extent  of  both  kinetic  and  chemical 
nonequilibrium,  the  propellant  gas  mixture,  and  the  arcjet  nozzle  geometry.  Other  factors  which 
to  the  overall  performance  include  the  thermal  loading  of  the  anode,2  and  the  voltage 
sneaths  at  the  electrodes.  A  comprehensive  model  of  an  arcjet  thruster  is  highly  complex, 

including  plasmadynamic,  fluid  dynamic,  radiative,  and  surface  and  volumetric  heat  transfer 
phenomena. 

Most  arcjet  design  strategies  to  date  have  been  empirical,  with  design  improvements  based 
on  experimental  observations.  The  need  for  a  more  complete  description  of  the  complex 
physical  processes  has  resulted  in  a  variety  of  numerical  efforts  ranging  from  simplified  one¬ 
dimensional  models  'o  to  more  comprehensive  two-dimensional  descriptions  of  the  gas 
dynamics  and  energy  transfer  processes.9"22  Comprehensive  arcjet  models  are  required  to  (1) 
interpret  experimental  results,  (2)  understand  the  physical  processes  in  regions  of  the  thruster 


where  diagnostic  techniques  are  challenging,  (3)  improve  arcjet  performance  and  (4)  scale  the 
devices  to  higher  or  lower  power  levels. 

II.  Nonequilibrium  Model 

An  axisymmetric,  steady,  laminar,  continuum  flow,  two-temperature  kinetic  nonequilibrium 
model  is  fomulated  for  a  direct  current  arcjet  with  flow  swirl,  a  variable  nozzle  geometry,  and  a 
variable  mixture  ratio  of  nitrogen  and  hydrogen.  Flow  swirl  is  modeled  by  including  an 
azimuthal  momentum  equation.  Two  methods  are  used  to  determine  species  populations.  If 
chemical  equilibrium  is  assumed,  the  species  populations  are  computed  from  a  generalized  law 
of  rnass  action.  For  chemical  nonequilibrium  simulations,  populations  are  determined  from 
species  continuity  equations. 

The  model  predicts  thiust,  specific  impulse,  and  internal  fields  for  pressure  p,  mass  density 
P’  seven  species  densities  ni,  and  anode  Ta,  electron  Te  and  heavy  species  Tg  temperatures. 
Additionally,  the  model  predicts  the  current  density  distribution  j,  voltage  potential  V,  and 
velocity  components  u,  v,  and  w. 

The  details  of  the  model  are  presented  in  six  sections  below  as;  (1)  the  fluid  dynamic 
equations,  (2)  the  energy  equations,  (3)  the  electromagnetic  equations,  (4)  the  species  number 
density  equations,  (5)  transport  coefficients,  and  (6)  boundary  conditions.  The  arcjet  geometry  is 
described  as  (x,r,6)  where  x  is  the  axial  coordinate,  r  is  the  radial  coordinate,  and  0  is  the 
azimuthal  coordinate. 

1.  Fluid  Dynamic  Equations 

The  viscous  fluid  dynamic  equations  are  summarized  below: 

Axial  momentum: 

0(puu)  ^  1  8(rpuv)  _4 

r  3  9^  9x/  <■  ^  drj  /i\ 

2  2  ,  I  dp 

3  3m  3r/  3r  3x  r  ^  3x/  3x 


Radial  momentum: 


3(puv)  1  3(rpvv)  _  ,  4. 

3x  r  9r  ^  9x/  3r  ^  drj 
+ +  3u  2v  gfi  4  M-v 
w  drj  3r  dx  3r  dr  3  j.2 

-  ^  I 

3r  ^  3x/  3r  r 


Azimuthal  momentum  (dp/dG  =  0): 


^M!I)).Li((^pv.2^Xw^))  = 

3x  r  9r 


(3) 


Axisymmetric  continuity: 


2.  Energy  Equations 

The  assumption  of  kinetic  nonequilibriuni  gives  separate  energy  equations  for  the  electrons 
and  heavy  species.  For  the  chemical  equilibrium  model,  the  electron  energy  equation  is: 

V.((neV-i)(CpeTe))  =  V.(AeVTe) 

Cf  /  ^ 

+  ^  Ij^^CpeTeDeVyej  +  -  elst  -  radiation 

where  'elst'  represents  energy  lost  through  elastic  collisions  with  the  heavy  a)ecies,23  and 
radiation  denotes  optically  thin  radiation  loss  due  to  continuum  bremsstrahlung.2^  The  elastic 
energy  transfer  source  term  is  calculated  as  the  sum  of  electron-molecule,  -atom,  and  -ion 
contributions  as: 


elst  =  3kb(Te-Tg)neme  I  ^  (6) 

i  e  "’i 

where  the  average  collision  frequencies  v^j  between  electrons  and  heavy  species  are  calculated 
using  the  mean  electron  thermal  speed  and  collision  cross  sections  Qes  as:^^ 

O) 

where  the  electron-ion  cross  sections  are: 


-  ice^kc  InAe 
2(kbTe)2 


(8) 


Elytron-molecule  collision  frequencies  are  multiplied  by  an  energy  loss  factor  6  to  account  for 
inelastic  losses  to  the  internal  energy  modes.  8.  25  a  constant  loss  factor  of  5  =  3000,  chosen  for 
tne  simulations  presented  here,  is  discussed  below. 

The  chemical  equilibrium  energy  equation  for  heavy  species  (molecules,  atoms,  and  ions)  is: 


where 


V  {nhvH)  =  V  (A-hVTg) + 

I  ^v  (p(CpiTg  +  hi°)  Di  Vyi )  +  elst  +  V  {v  •  x  ) 


(9) 


(ao) 

The  ionization  and  dissociation  energies  of  molecules  and  atoms  are  arbitrarily  included  in  the 
nea^  species  energy  eqimtion  as  reference  enthalpies  hi^.  The  reference  enthalpies  and  specific 
neats  are  summanzed  in  Table  1.  t'  t' 

For  chemical  nonequilibrium  simulations,  a  finite  rate  chemistry  model  is  used  to  determine 
inelastic  energy  exch^ge  rates  between  the  particle  classes.  This  model  is  summarized  in  Table 
z,  wnicn  lists  the  individual  reaction  rates  and  energy  exchanges  assumed  for  each  chemical 
process.  1  he  species  continuity  equations  are  multiplied  by  the  reference  enthalpies 
associated  with  each  chemical  reaction,  and  the  terms  appearing  in  the  energy  equations  as  the 
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convection  and  diffusion  of  ionization  and  dissociation  potentials  are  replaced  with  energy 
transfer  rates.  The  electron  energy  equation  for  chemical  nonequilibrium  is  then 

V.((neV-J-)(CpeTe))  =  V.(XeVTe) 

+  ^  Ij^^CpeTeDeVyej  +  -  elst  -  radiation  -  X  Le 

where  individual  energy  loss  rates  Le  are  summarized  in  Table  2. 

The  chemical  nonequilibrium  heavy  species  energy  equation  is: 

V  (nhvH)=V  (>.hVTg)+  I  -LV-(pCpiTg  DjVyi ) 

+  elst+  V  (v  •  T  )  -  XLh 

where  the  total  energy  is  now: 


The  anode  ternperature  distribution  is  determined  from  the  energy  equation  for  heat 
conduction  in  a  solid.  Neglecting  ohmic  dissipation  in  the  anode  gives 

V(XaVTa)  =  0  (14) 


(11) 


(12) 


where  =  /laCTa).^® 


3.  Species  Equations 

For  chemical  equilibrium  (CE)  simulations,  the  species  number  densities  are  given  as  a 
junction  of  the  pressure  p,  the  heavy  species  temperature  Tg,  and  the  electron  temperature  Te-  A 
^ven-species  plasina  composition  of  N2,  H2,  N,  H,  N"*",  H"*",  and  electrons  is  assumed. 
Dissociation  and  ionization  equations  are  constructed  from  a  generalized  law  of  mass  action  as 
given  by  entropy  maximization. 31  This  yields  the  following  equations  for  the  nitrogen  and 
hydrogen  species: 


(15a) 


(15b) 


(15c) 
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nenN+ 

"N 


(15d) 


In  the  above  expressions,  the  excitation  partition  functions  Zex  are  functions  of  Tg,  while 
Zrot  and  Zvib.are  functions  of  Tg.  Details  for  the  calculation  of  the  various  partition  functions 
can  be  found  in  References  24  and  3 1  -34. 

The  plasma  equation  of  state,  quasi-neutrality,  and  the  definition  of  the  mixture  parameter 
Xmol  provide  the  additional  relations  required  to  solve  for  the  seven  unknown  species.  Note  that 
the  initial  mixture  composition  is  written  as  XmolN2  +  H2  so  that  pure  hydrogen  (xmol  =0), 
simulated  hydrazine  (xmol  =  1/2)  and  simulated  ammonia  (xmol  =  1/3)  can  be  easily 
investigated:  ^ 


_  2nN2  +  HN  +nisj+ 

^moi  — - - 

2nH2  +  "H  +nH+ 


(16) 


ne  =  nH+  +  nN+  (17) 

p  =  nekbTg  +  (nNj  +  nH2  +  iN 
+  nH  +  ne)kbTg 

For  chemical  nonequilibrium  (CNE)  simulations,  Eqs.  (15)  and  (16)  are  replaced  by 
individual  species  continuity  equations,  which  are  written  in  terms  of  the  species  mass  fractions 
yi  =  mini/p: 


V(pvyi)  =  V(pDiVyi)+  p; 


(19) 


..Mixture-averaged  in^s  diffusion  coefficients  Di  are  employed,  rather  than  multi-component 
diffusion  coefficients.  With  this  method,  Eq.  (19)  is  solved  for  yN2>  YHj  YN.  YN+j  and  yH+-  TTie 
electron  population  is  then  calculated  from  Eq.  (17),  and  yH2  is  determined  from  overall  mass 
continuity,  or  equivalently  Zyi  =  1.  The  species  productions  rates,  pi  =  mini,  are  determined  as  a 
runcuon  of  Tg  and  Tg  from  the  finite-rate  chemistry  processes  summarized  in  Table  2.  The 
reterences  tor  the  individual  reaction  rates  are  summanzed  in  the  right  column  of  Table  2.  For 
two-™  reactions,  the  reverse  reaction  rates  are  extracted  from  the  forward  reaction  rate  and  the 
equilibnum  constants'  given  by  Eqs.  (15a-d). 


4.  Electromaenetic  Equations 

The  calculation  of  the  electrical  current  distribution  is  critical  to  understanding  arc  physics 
and  poses  a  chafiengmg  problem  in  arcjet  analysis.  The  current  distribution  is  largely 
dete^ned  by  the  flow-field  electrical  conductivity  a,  which  is  a  strong  function  of  the  electron 
nunaber  density  and  temperature  distribution.  Models  which  assume  local  thermodynamic 
equilibrium  (LTE)  must  artificially  elevate  a  in  the  cool  boundary  layer  region  to  avoid 
excessive  voltage  drops  near  the  anode.  11.15  xhe  artificial  restriction  that  Tg  =  Tg  results  in  an 
underprediction  m  the  ionization  (and  therefore  electron  number  density  and  o)  in  the  boundary 
myer  regions  and  a  subsequent  overprediction  in  power  deposition.  Hence,  the  'conductivity 
^^,?^PlOY®d  elsewhere iii’i'+  to  provide  agreement  between  the  experimentally-observed 
Sdd  predicted  power  depositions. 

Diffusion  of  electrons  from  the  arc  core  towards  the  anode  wall  also  elevates  the  electron 
nurnber  density  and  plays  an  important  role  in  determining  the  boundary  layer  electrical 
conductivity.  This  has  been  demonstrated  by  Butler,  Kull,  and  King,  1 1  who  show  that  the 

diffusion  provides  a  reasonable  current  distribution  without  imposing 
artiricial  restrictions  on  cr  or  the  current  attachment  region.  However,  the  two-temperature 
solutions  of  Miller  and  Martinez-SanchezI'  and  Keefer  et  al.^6  demonstrate  that  significant 
^netic  nonequihbnum,  with  elevated  Tg,  exists  outside  of  the  arc  core  in  the  constrictor  region, 
inese  results  indicate  that  both  diffusion  and  increased  ionization  due  to  elevated  Tg  are 
important  factors  which  determine  the  boundary  layer  electrical  conductivity. 
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The  correct  approach  to  the  arc  attachment  problem  is  to  model  the  anode  as  an  equipotential 
surface  and  allow  the  current  distribution  to  be  a  model  output  which  is  independent  of  artificial 
restrictions.  This  is  a  self-consistent  approach  to  solving  for  the  current  distribution  and  electric 
field  simultaneously,  for  which  the  current  distribution  is  properly  coupled  to  the  flow-field. 

Assuming  induced  magnetic  fields  are  negligible,  and  neglecting  electron  pressure  gradient 
driving  terms.  Ohm's  law  and  current  conservation  are  used  to  solve  for  j  and  V:25 


.VV  =  j/a 

(20) 

v-j  =  0 

(21) 

5.  Transport  Coefficients 

Transport  coefficients  for  this  model  are  calculated  using  mean  free  path  mixture 
rules.'^^  The  viscosity  is  computed  as: 


11=  I  — lii  (22) 

I  njMij 
j=l 


where  Mjj  is  defined  as 


Mij  = 


- A<2’2) 

2niij  “ij 


mi  -K2,2) 


(23) 


and  the  pure  species  viscosity  is 


=  l/2miniCili 


(24) 


where  ci  is  the  mean  thermal  speed,  li  is  the  mean  free  path,  and  is  the  energy- 
averaged  collision  integral  for  interaction  between  species  i  and  j. 

Ths  ^anslational  contribution  to  the  heavy  species  thermal  conductivity  is 
calculated  using  an  equation  similar  to  (24)  with  the  pure  species  conductivity  given  as: 


9Yi-5 
ISyj  - 15 


(25) 


The  electron  thermal  conductivity  is: 


I  njMjj 

j=i 


and  electrical  conductivity  is  calculated  from: 


(26) 
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The  effective  diffusion  coefficients  Di  are  calculated  from: 


Di  =  — 

I  xj/Dij 
J^i 

where  the  Dij  are  the  binary  mixture  diffusion  coefficients  given  as:^^ 


3 

pkbT(mi+  mj) 

loniljj 

\  TCrnjmj  j 

(28) 


(29) 


The  electron  and  ion  diffusion  coefficients  are  modified  to  account  for  the  ambipolar 
electric  field.  Assuming  =  Z!1h+ =  ZI1N+ ,  the  ambipolar  diffusion  coefficients  for 

"e  nH+  nN+ 

electrons,  H+,  and  N+  are: 


^^e(M-H+"H+  +  I^N+"N+) 


+  (DH+nH+ +  J 

M-e"e  +  ^^H+nH+  +  ^^N+"N+ 

(30) 

*^H+(^^e"e  +  I^N+”N+) 

+  ^^H+  (DeHe  "  DN+nN+) 

(31) 

M-e"e  +  ^^H+”H+  +  1^N+"N+ 

I^N+(^^e"e  +  M^H+"H+) 

+  M^N+  (DeHe  -  Dn+nm.) 

Me  +  f^H+"H+  +  M+"N+ 

(32) 

where  the  mobilities  p  are  calculated  from 

Ti 

.  The  above  relations  require  the  energy-averaged  collision  integrals  for  28  species 
interactions.  Collision  integrals  and  collision  cross  sections  are  obtained  from 
Keterences  37-43.  Temperature-dependent  curve  fits  are  used  for  the  transport 
property  and  elastic  transfer  calculations. 

6.  Boundary  Conditions 

The  physical  boundary  conditions  for  the  continuity,  momentum,  energy,  and 
electromagrietic  equations  are  summarized  in  Table  3  and  Fig.  1.  At  the  inflow 
Dounda^  the  flow  is  axial  and  the  total  temperature  Tt  and  total  pressure  pt  are 
specified.  Tf^  total  pressure  is  adjusted  to  obtain  the  desired  mass  flow  rate.  Ine  u- 
veiocities  and  the  inlet  static  temperatures  are  then  computed  from  the  isentropic 
stagnation  property  relations.  The  effects  of  flow  swirl  are  mcorporated  by  specifying 
the  azimuthal  mlet  w-velocities  as  a  fraction  of  the  u-velocities.  At  the  solid  cathode  1C 
and  anode  ID  surfaces,  no-slip  conditions  are  employed.  Zero  radial  gradients  are 
specified  on  the  center-line  IE. 


4.U  (^tflow  boundary  IB,  the  static  pressure  is  extrapolated  from  the  interior  of 

the  flow.  The  outflow  U;-velocities  are  updated  to  satisfy  continuity  using  velocity 
corrections  as  given  by  discretized  forms  of  the  momentum  and  continuity  equations 
for  the  exit  plane.  This  is  consistent  with  the  numerical  solution  method  discussed  in 
section  III  below.  Note  that  a  zero  electron  temperature  gradient  is  employed  at  the 
walls  1C  and  ID.  At  low  temperatures,  the  electron  thermal  conductivity  approaches 
zer^o  that  this  boundaiy  condition  will  have  little  effect  on  the  results.44 

The  boundary  conditions  for  the  voltage  are  shown  in  Table  3.  The  anode  is 
modeled  as  an  equipotential  surface,  while  the  cathode  voltage  is  equipotential  with  the 
voltage  updated  to  obtain  the  specified  total  current  I.  Note  Qiat  in  this  formulation  the 
electrode  sheath  voltage  drops  Vs  have  been  uncoupled  from  the  bulk  plasma.  The 
axial  gradient  of  the  voltage  is  set  to  zero  at  the  inflow  and  outflow  boimdaries  so  that 
current  is  forced  to  attach  within  the  arciet  nozzle.  Recent  radial  electric  field  and  <7 
measurements  along  the  exit  plane,  showing  the  existence  of  roughly  1  A/cm^  current 
den^ty,  suggest  that  this  boundary  condition  may  need  refinement.^S 

The  boundary  conditions  for  the  anode  energy  equation  are  also  shown  in  Table  3. 
On  the  outer  surfaces  2B  and  2D,  the  local  heat  flux  is  given  ^  radiation  to  the 
surroundings  at  a  specified  background  temperature  Tbe  <<  Ta  and  a  constant 
emissiyity  of  e  =  0.31.  Along  the  anode  surface  2c:,  the  net  heat  flux  is  calculated  as  the 
sum  of  (1)  the  heat  flux  due  to  conduction  and  species  diffusion  from  the  flow-field,  (2) 
the  radiation  loss  through  the  exit  plane  to  the  surroundings,  and  (3)  the  sheath  losses 
I’Y  the  product  oTthe  local  current  density  and  an  assumed  sheath  voltage  drop 
•  r  •  Note  that  radiation  exchange  between  the  plasma  and  the  anode  is  neglected 
in  mis  lormulation.  The  upstream  anode  surface  2A  is  either  a  fixed  temperamre  or 
zero  heat  flux  boundary.  ^ 

The  chemical  nonequilibrium  model  requires  boundary  conditions  for  the  species 
populations.  At  the  anode  surface,  the  diffusion  velocities  are  equated  to  a  fracfion  ai 
of  the  species  thermal  velocities:29  ^ 


3n  4 


Kmi 


(33) 


pie  anode  is  assumed  to  be  catalytic  for  recombination,  so  that  aH+  =  aN+  =  aH  =  ocm  = 
1.  Assuming  that  molecules  are  reflected  from  the  surface  gives  aul  =  «N2  =  0.  Ide 
inflow  IS  in  chemical  equilibrium  up  to  the  constrictor  entrance,  where  the  chemical 
nonequilibrium  srmulaiion  begins.  The  mass  fractions  at  the  exit  plane  are  extrapolated 
from  me  mtenor  of  the  flow.  ^ 


III.  Numerical  Solution  Method 


The  solu^ton  method  is  structured  around  the  global  continuity  and  momentum 
equations,  .pese  are  solved  using  a  compressible  form  of  the  pressure-based  PISO 
algorithm,  where  the  density  variatioris  are  implicitly  included  in  the  pressure- 
correction  procedure  as  described  by  Rhie.^'  Staggered  velocity  and  scalar  grids  yield 
a  well-coimected  pressure  field.  Upstream  weighted  densities  are  employed  to  insure 
numerical  stability.'*  The  equations  are  discretized  over  finite  cell  volumes,  and 
power  law  differencing  of  convective  and  diffusive  portions  of  the  linearized 
coefficients  Reserves  numerical  stability  in  regions  of  the  flow  where  convection 
dominates.  The  solution  process  is  iterative,  and  values  for  the  field  variables  from  the 
latest  iteration  are  used  to  linearize  the  source  terms.  Under-relaxation  is  employed  so 
solutions  are  fractionally  added  to  the  solution  from  the  previous 
Iteration.  1  he  govemmg  equations  are  transformed  into  natural  coordinates  and  solved 
on  a  imiform  computational  mesh. 

Figure  2  illustrates  the  solution  algorithm.  The  energy  equations  are  solved  to 
update  gas,  electron,  and  anode  temperature  fields.  The  species  number  densities, 
transport  coefficients,  and  source  terms  are  then  updated  using  the  latest  available 
temperatures  and  pressures.  The  continuity  and  momentum  equations  are  solved  to 
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update  the  pressure  and  velocity  fields.  Finally,  the  voltage  and  current  density 
distributions  are  updated.  Iterations  are  continued  until  (1)  the  normalized  energy 
equation  residuals  are  within  1%  of  the  electrical  power  input,  (2)  the  difference 
between  the  inlet  and  exit  plane  mass  flow  rates  is  less  than  1%  of  the  total  mass  flow 
rate,  and  (3)  the  axial  and  radial  momentum  equation  residuals  are  within  1%  of  their 
respective  exit  momentum  fluxes. 

IV.  Results  and  Discussion 


Overview  of  Previous  Simulations 

Before  proceeding  to  the  present  results,  a  discussion  of  the  effect  of  the  electron- 
X?  A  c  A^T  ^oss  factor  5  on  our  initial  voltage  predictions  is  appropriate.  For  the 

^  kW  arcjet  geometry,  we  found  that  the  6-factors  quoted  in  Ref.  25  (6  == 
U(lU-iUU)j  resulted  m  current  attachment  in  the  constrictor,  with  corresponding 
underpredictions  of  the  operating  vqltage  of  V  =  40  V  at  I  =  10  A.  We  also  noted  that 
increasing  the  mass  flow  rate  had  little  or  no  effect  on  the  attachment  location.  We 
electrically  msulated  the  constrictor  and  observed  little  effect  on  the  voltage  prediction 
Martinez-Sanchez  recently  re-examined  the  6-factors  for  N?,  and  values  as  high  as  6  = 
** -iu  noted.®  We  choose  6  =  3000  because  (1)  reasonable  convective  effects  result, 
wim  the  arc  bemg  swept  further  downstream  as  the  mass  flow  rate  is  increased,  as  was 
oriemally  postulated  in  Ref.  8,  and  (2)  better  agreement  is  achieved  between  the 
predicted  and  experimentally  observed  operating  voltage,  while  maintaining  the 
equ^otential  anode  voltage  boundary  condition  in  the  constrictor. 

The  constrictor  has  been  electrically  insulated  ^  others  to  provide  better  agreement 

predicted  voltage. These  models  simulated  the  German 
^  ^  j  i  ^  higher  power  10  kW  arcjet  with  an  =  5  mm  long  constrictor,  and 

predicted  that  nearly  all  of  the  electrical  power  was  deposited  in  the  constrictor,  l^e 
predictions  for  our  geometry,  where  Xcon  =  0.25  mm,  indicate  that  roughly  50%  of  the 
power  IS  deposited  downstream  of  the  constrictor. 

Results 

We  present  Aree  simulations  for  the  geometry  of  a  NASA-Lewis  1  kW  arcjet.  These 
are  summarized  m  T^ble  4,  where  mass  now  rate,  voltage,  current,  thrust,  and  specific 
impulse  are  shown.  The  simulations  are  designated  as  either  chemical  equilibrium  (CE) 
or  chemical  nonequilibrium  (CNE).  For  each  case,  the  upstream  anode  surface  2A  (Fig. 

ii«  7  1000  K,  and  the  cathode  1C  temperature  varies 

hnearly  from  Tc  —  1000  K  at  the  inlet  to  Tc  =  3000  K  at  the  cathode  tip.  The  constrictor 
diameter  is  0.63  nam.  The  cathode  gap  spacing,  as  measured  axially  from  the 
convergmg  portion  of  the  anode  surface,  is  Xcat  =  0.58  mm  for  the  hydrazine  cases  and 
the  hydrogen  case.^M  anocfe  voltage  drop  of  Vo  =  13.0  V  is  assumed, 
aim  the  anode  emissivity  is  e  =  0.31.12  jj^e  inlet  flow  swirlvelocity  is  30%  of  the  axial 
imiow  velocity,  pie  computational  grids,  shown  in  Fig.  3,  are  comprised  of  62  axial  by 
1/  radial  nodes  for  the  plasma,  and  62  axial  by  9  radial  nodes  for  the  anode  domain, 
bimuiation  run  time  is  several  hours  on  a  Convex  C-240  mainframe  computer. 

1.  Hydrazine  Simulations 

pie  CE  and  CNE  hydrazine  simulations  are  systematically  compared  in  Figs.  4  - 15. 
^ontours  ot  cunent  density  and  anode  temperature  are  indicated  in  Fig.  4  for  CE 
Hydrazine.  At  the  cathode  tip,  the  current  density  is  j  -  40,000  A/cm2.  A  peak  anode 
rrt  density  of  j  =  18  A/cm2  is  noted  at  a  location  of  x  «  11  mm,  or 

ougmy  5  mrn  downstream  of  me  constrictor.  The  anode  temperature  increases  from 
boundary  value  of  Ta  =  1000  K  to  Ta  «  1300  K  near  the  exit  plane.  The 
'rNt  Iwdrazme  results  are  shown  in  Fig.  5.  The  peak  anode  current  density  of  j  «  15 
/cm  occurs  at  x  ~  8  mm,  upstream  from  that  for  the  CE  case.  The  anode  temperature 
is  somewhat  lower  than  for  CE,  with  a  value  of  Ta  *  1200  K  at  the  exit  plane. 

•1,  ^^^ked  difference  in  the  current  attachment  predictions  is  more  clearly 
Illustrated  in  Fig.  6,  where  the  anode  current  densities  are  indicated.  For  CNE,  a  more 
dittuse,  bimodm  distribution  is  noted,  with  the  attachment  closer  to  the  constrictor, 
^consequently  the  voltage  prediction,  indicated  in  Table  4,  is  significantly  lower  at  87  V 
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for  CNE  versus  113  V  for  the  CE  simulation.  The  anode  Tg  distributions  are  indicated 
^  nonequilibrium  is  noted,  with  maximum  electron 

re^  ectiv^^^^^  ~  ^  ~  ^  for  the  CE  and  CNE  simulations 

Heavy  species  and  electron  temperature  contours  in  the  constrictor  region  are 
mdicated  mEi^  8-11  In  Fig.  8,  a  maximum  Tg  =  22,000  K  occurs  near  the  center  of  the 
constrictor  for  CE  hydrazine.  The  central  regioh  of  the  arc  is  near  thermal  equilibrium 
as  mdicated  by  a  maximum  Tg  =  23,500  K  shown  in  Fig.  9.  The  result  for  CNE 
FigurS^O  and  maxima  of  Tg  «  20,000  K  ancTTg  =  23,000  K  indicated  in 

yeloci^  contours  are  shown  in  Figs.  12  and  13.  A  maximum  centerline 
velocity  ot  u  »  8  xm/ s  occurs  at  x  =  8  mm,  and  then  decreases  because  of  heating  and 
viscpus  losses  to  u  ~  6  km/s  at  the  exit  for  CE  hydrazine.  For  CNE  hydrazine,  a 
maximum  velocity  of  u  «  7  km/ s  is  predicted  to  occur  at  x  =  9  mm,  decreasing  to  u  «  6  5 
km /s  at  the  exit.  ° 

Finite-rate  chemistry  must  be  included  to  capture  frozen  flow  effects.  In  Figs.  14 
and  15,  exit  plane  predictions  of  electron  number  density  ng  are  presented.”  The 
predicted  centerlme  density  for  CE  hydrazine  is  ng  =  1.6  x  10^1  cm-3,  while  for  CNE,  the 
density  IS  3  orders  of  magnitude  higher  at  ng  =  1.6  x  10^^  cm*^.  Also  shown  are  Tg  and 
thermal  nonequifibrium  is  indicated  for  the  CE  case  with  Tp/To- »  2.0 
arthe  centerlme,  while  the  CNE  case  is  near  equilibrium  at  Tg/Tg  =  1.1.  ° 

2.  l^drogen  Simulation 

The  results  for  the  CNE  hydrogen  simulation  are  indicated  in  Figs.  16  -  20.  Anode 

density  contours  are  shown  in  Fig.  16.  The  peak  anode 
current  density  of  j  =  16  A/cm^  occurs  at  the  attachment  location  of  x  =  8  mm.  Due  to 
tne  nigner  specific  power  deposition,  the  thermal  loading  of  the  anode  is  somewhat 
1  ^^^l\ydxazine  cases,  with  a  predicted  temperature  near  the  exit 

plane  of  Ta  ~  1400  K.  Profiles  of  Tg  and  j  at  the  anode  surface  are  shown  in  Fig.  17.  The 
results  are  qualitatively  similar  to  the  CNE  hydrazine  case,  with  high  thermal 
nor^quihbriuin  of  Tg  =  12,000  K  predicted  in  the  current  attachment  region.  ” 
r-  1  species  and  electron  temperature  contours  are  shown  in  Figs.  18  and  19.  In 
Fig.  18,  a  peak  temperature  of  Tg  =  13,500  K  is  noted.  The  electron  temperature  is 
similar  in  the  central  region  of  thl  flow  with  Tg  =  14,000  K  indicated  in  Fig.  19.  At  a 
^700  ^  n  C  Tg  decreases  in  the  radial  direction  reaching  a  minimum  of  Tg  = 

£  T  ~  ^jm.  The  electron  temperature  then  increases,  eventually  reaching  a 
^  ~  12,000  K  at  the  anode  surface.  The  electron  temperature  profile  is 

qualitatively  similar  to  the  nonequilibrium  predictions  of  Miller  and  Martinez- 
banchez^'  and  Keefer  et  al.i” 


Axial  velocity  contours  are  shown  in  Fig.  20.  A  maximum  centerli 
km/s  occurs  at  x  —  7  mm  decreasing  to  u  =  12  km/s  at  the  exit  plane. 


centerline  value  of  u  »  15 


V.  Conclusions 


detailed,  propellant-flexible,  N2/H2  arqet  model  is  developed,  and  is  sufficiently 
used  for  arcjet  design  studies.  A  thermal  nonequilibriunu 
cnemical  equilibrium  model  was  first  dev^oped,  and  then  generalized  to  chemical 
nonequilibrium  by  including  finite  rate  chemistry  and  mass  diffusion.  The 


modifications  to  the  electrical  conductivity. 

1  y^®^^^^|^®builibrium  and  nonequilibrium  results  have  been  compared  for 
J^arazme.  High  :^ar-anqde  thermal  nonequilibrium  is  indicated  in  both  cases  with 
Af/  T  =  12  for  CE  and  CNE  predictions  respectively.  In  the  constrictor, 

me  central  region  ofthe  arc  is  near  thermal  equilibrium  with  Tg  =Te  =  20,000  K.  Mass 
dittusion  and  finite  rate  chemistry  dramatically  affect  the  ale  structure.  The 
incorporation  of  these  chemical  nonequilibrium  processes  reduces  the  arc  length  and 
lowers  the  predicted  operating  voltage.  Additionally,  higher  frozen  flow  losses  are 
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npted.  Predicted  electron  densities  at  the  ardet  exit  plane  are  three  orders  of  magnitude 
higher  with  nonequilibrium  chemistry.  ° 

Chemical  npnequilibrium  results  have  been  presented  for  hydrogen.  The  results  are 
qualitatiyeiy  similar  to  those  for  Iwdrazine.  However,  lower  arc  temperatures  are 
with  interline  Te  »  Tg  »  14,000  K  in  the  constrictor. 

The  combined  effects  of  viscous  losses  and  supersonic  heating  were  noted  in  all 
cases.  Centerlme  axial  velocities  reach  a  maximum  a  few  millimeters  downstream  of 
the  constrictor,  and  then  decrease  to  exit  plane  values  of  u  ==  6  km/s  and  u  =  12  km/s 
tor  hydrazme  and  hydrogen  respectively. 

An  anode  temperature  model  was  mcluded.  For  CE  hydrazine,  maximum  anode 
temperatures  of  Ta  *  1300  K  are  predicted  near  the  exit  plane.  For  CNE  hydrazine  and 
maxima  of  Ta  »  1200  and  Ta  «  1400  are  predicted. 

The  model  has  generated  much  useful  data;  however,  vmresolved  issues  remain. 
1  he  voltage  predictions  largely  depend  on  the  arc  attachment  location.  We  have  found 
that  parametncally  mcreasmg  electron-molecule  inelastic  energy  transfer  increases  the 
arc  length.  A  more  general  physical  model  of  electron-molecule  inelastic  transfer  is 


predict  i 


/oltage  drops. 


completG  description  of  the  coupling  between  me  electron  energy  and  molecular 
internal  energy  modes.  This  maybe  required  to  describe  the  arc  attachment  and  near¬ 
anode  electrical  conductivity  more  accurately. 
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Table  1:  Species  reference  enthalpies  and  specific  heats. 
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Figure  1:  Schematic  of  arq'et  thruster  indicating  fluid  and  anode  domain  boundaries. 
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Table  3:  Boundary  conditions  for  the  fluid  dynamic,  energy,  and  electromagnetic 
equations.  Energy  equation  boundary  conditions  are  shown  for  both  the  fluid  and  anode 
domains.  For  radiative  boundary  conditions,  Ts  is  the  local  anode  surface  temperature. 
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"^ble  4:  Summary  of  predictions.  CE  denotes  chemical  equilibrium  and  CNE  denotes 
chemical  non-equilibrium  simulations. 


39 


are  residuals  for  all 
conservation  equations 
converged? 


yes 


stop 


start 


initialize  thruster  geometry 
field  variables,  transport 
coefficients,  energy  sources 


Figure  2:  Numerical  solution  algorithm.  The  solution  method  is  iterative,  employing 
successive  imder-relaxation  of  the  indicated  field  variables. 


cnnstrictor  region 


Figure  3:  Computational  grids.  The  grid  is  tuned  axially  for  high  density  in  the 
constrictor  and  radially  for  high  density  near  the  electrode  surfaces. 
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Figure  4:  Anode  temperature  and  current  density  for  hydrazine;  chemical  equilibrium 
(CE)  simulation.  The  maximum  anode  current  density  is  j  =  18  A/cm^  at  x  ~  11  mm. 
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Figure  5:  Anode  temperature  and  current  density  for  hydrazine;  chemical 
nonequilibrium  (CNE)  simulation.  The  maximum  anode  current  density  is  j  ~  15 
A/cm2  atx»8mm. 
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Figure  6:  Current  density  at  the  anode  surface  for  hydrazine.  Results  are  shown  for 
chemical  equilibrium  (CE)  and  chemical  nonequilibrium  (CNE)  simulations. 
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Figure  7:  Electron  temperature  at  the  anode  surface  for  hydrazine.  Results  are  shown 
for  chemical  equilibrium  (CE)  and  chemical  nonequilibrium  (CNE)  simulations. 
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Figure  8:  Heavy  species  temperature  for  hydrazine;  chemical  equilibrium  (CE) 
simulation.  The  maximum  is  Tg  =  22,000  K  at  the  constrictor  centerline.  The  exit  plane 
is  located  at  x  =  18.5  mm. 
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Figure  9:  Electron  temperature  for  hydrazine;  chemical  equilibrium  (CE)  simulation. 

The  maximum  is  Te  =  23,500  K  at  the  constrictor  centerline.  The  exit  plane  is  located  at 
X  =  18.5  mm. 
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Figure  10:  Heavy  species  temperature  for  hydrazine;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  Tg  ~  20,000  K  at  the  constrictor  centerline.  The  exit  plane 
is  located  at  x  =  18.5  mm. 
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Figure  11:  Electron  temperature  for  hydrazine;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  Te  =  23,000  K  at  the  constrictor  centerline.  The  exit  plane 
is  located  at  x  =  18.5  mm. 
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Figure  12:  Axial  velocity  contours  for  hydrazine;  chemical  equilibrium  (CE)  simulation. 
The  maximum  is  u  *  8  km/ s  at  x  =  8  mm.  The  centerline  velocity  decreases  to  u  »  6 
km/ s  at  the  exit  plane. 
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Figure  13:  Axial  velocity  contours  for  hydrazine;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  u  «  7  km/ s  at  x  =  9  mm.  The  centerline  velocity  decreases 
to  u  «  6.3  km/ s  at  the  exit  plane. 
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Figure  14:  Radial  distributions  of  electron  density,  electron  temperature,  and  heavy 
species  temperature  at  the  exit  plane  for  hydrazine;  chemical  equilibrium  (CE) 
simulation. 
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Figure  15:  Radial  distributions  of  electron  density,  electron  temperature,  and  heavy 
species  temperature  at  the  exit  plane  for  hydrazine;  chemical  nonequilibrium  (CNE) 
simulation. 
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Figure  16:  Anode  temperature  and  current  density  for  hydrogen;  chemical 
nonequilibrium  (CNE)  simulation.  The  maximum  anode  current  density  is  j  =  16 
A/ cm2  at  X  «  8  mm. 


Figure  17:  Electron  temperature  and  current  density  at  the  anode  surface  for  hydrogen; 
chemical  nonequilibriiun  (CNE)  simulation. 
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Figure  18:  Heavy  species  temperature  for  hydrogen;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  Tg  =  13,500  K  at  the  centerline  near  the  cathode  tip.  The 
exit  plane  is  located  at  x  =  18.5  mm. 
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Figure  19:  Electron  temperature  for  hydrogen;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  Te  =  14,000  K  at  the  centerline  near  the  cathode  tip.  High 
near-anode  thermal  nonequilibrium  is  indicated  downstream  of  the  constrictor,  with  Te 
=  12,000  K  at  the  anode  surface.  The  exit  plane  is  located  at  x  =  18.5  mm. 
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Figure  20:  Axial  velocity  contours  for  hydrogen;  chemical  nonequilibrium  (CNE) 
simulation.  The  maximum  is  u  »  15  km/s  at  x  «  8  mm.  The  centerline  velocity 
decreases  to  u  «  12  km/ s  at  the  exit  plane. 
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